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Recent developments in the spectropolarimetric study of compact objects, specifically black 
holes (stellar and massive) and neutron stars are reviewed. The lectures are organized around 
five topics: disks, jets, outflows, neutron stars and black holes. They emphasize physical mech- 
anisms and are intended to bridge the gap between the fundamentals of polarimetry and the 
phenomenology of observed cosmic sources of polarized radiation, as covered by the other lec- 
turers. There has been considerable recent progress in spectropolarimetry from radio through 
optical frequencies and this is producing some unique diagnostics of the physical conditions 
around compact objects. It is argued that there is a great need to develop a correspondingly 
^ ' sensitive polarimetric capability at ultraviolet through 7-ray energies. 
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Spectropolarimetric observations, particularly those at radio and optical wavelengths, 



have played an important role in high energy astrophysics. From the discovery of syn- 
chrotron radiation to the first good evidence for AGN unification, from the polarization 
patterns in the coherent emission of radio pulsars to the discovery of variable, linear po- 
larization in the absorption troughs of broad absorption line quasars, polarization studies 
often provide the best and sometimes the only clue we have as to the geometric dispo- 
sition of the emitting elements in these diverse sources when we cannot resolve them 
directly. 

These notes summarize lectures delivered by Roger Blandford at the XII Canary Is- 
lands Winter School on Astrophysical Spectropolarimetry. They are written up with 
the assistance of Eric Agol (Disks), Leon Koopmans (Jets), Hee-Won Lee (Outflows), 
Jeremy Heyl (Neutron Stars) and Avery Broderick (Black Holes) The lectures were in- 
tended to provide a bridge between the general physical foundations of polarimetry and 
its practical description presented at the school by Drs. Landi DcgPInnocenti and Keller 
and the observationally oriented lectures of Drs. Antonucci and Hildebrand. They also 
make some important connections to solar, stellar (especially white dwarf) and maser 
polarimetry as described by Drs. Stenflo, Mathys and Elitzur, respectively. They are 
organised around five generic astrophysical sources: disks, jets, winds, neutron stars, 
and black holes. In each case a cursory motivation is provided by summarizing some 
relevant observations and presenting some of the key issues that polarimetry can help to 
resolve. This is followed by an heuristic discussion of some relevant physical mechanisms 
in a manner which, it is hoped, will allow them to be applied elsewhere followed by a 
brief account of how they have been used so far and some suggestions of possible future 
investigations. 



1. Disks 

1.1. Motivation 

Accretion disks are commonly found when gas, with angular momentum, is gravita- 
tionally attracted towards a central massive body (Frank, King, & Raine 1992). First 
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described in the context of Laplace's nebular hypothesis and first seriously analyzed by 
Lust (1952) they have been observed around black holes, neutron star and white dwarf 
binaries (Shapiro & Teukolsky 1983), around protostars, and especially within active 
galactic nuclei, including quasars (Blandford, Netzer, & Woltjer 1990, Krolik 1999). It 
is this last type of disk that provides us with much of our most detailed observations. 
Although it is not part of my task to discuss them, observations of young stellar objects 
and cataclysmic variables are turning out to be particularly instructive and much of what 
follows is informed by the results of these studies. 

Disks are planar structures and if their opacity is predominantly scattering, by cither 
free electrons or dust grains, then the direction and strength of the polarization tells us 
about the orientation and inclination of the disk as well as the location of the continuum 
source. As we shall discuss, (cf also Antonucci, Hildebrand, these proceedings), most 
astrophysical disks are associated with jets or bipolar outflows and when these can be 
resolved, they may represent the projected rotation axis of the inner disk. (As we shall 
also see, disks can be warped and this axis can change with radius and, in the case of 
AGN, it may be quite different from the axis of the host galaxy.) 

We wish to use polarization observations to determine what disks are really like. Un- 
fortunately, the current observational capability is limited. Polarimetry in the radio, the 
near infrared and the optical regions of the spectrum is really quite good by astronom- 
ical, (although not solar), standards. Optically, spectropolarimetry has been performed 
at the 0.001 level down to R ~ 18. Measurement in the mid and far infrared is more 
of a challenge, though one that has been met in the far infrared, (Hildebrand, these 
proceedings) . 

However, to understand the inner disk we need ultraviolet and X-ray polarimetry. The 
former was carried out for a while on bright quasars using the HST Faint Object Camera, 
as we shall describe in section 3 below, although this has proved to be a little controversial. 
X-ray polarimetry has really only been accomplished successfully on a few bright sources 
(Meszaros et al. 1988). There are plans to fly a more sensitive polarimetcr on Spectrum-X. 
It will become clear that there is a very strong scientific case to be made to develop X-ray 
polarimetry There is also a strong incentive to develop a 7-ray polarimetric capability 
though, here, the technical challenges are even greater. In principle, Compton telescopes 
operating at ~ McV energies, record polarimetric information though, in practice, it has 
proven to be almost impossible to extract this signal from observations taken to date. 

1.2. Observation 

A particularly good example of an AGN accretion disk can be found in NGC 4258 
(cf Elitzur, these proceedings). Here water maser observations reveal a resolved, disk 
orbiting a 43 million solar mass black hole (beyond all reasonable doubt). We have 
believed for a long while that gas moves radially inwards through this disk as a result 
of a hydromagnetic torque that transports angular momentum outward. The binding 
energy that is released by the infalling gas can be radiated away and this process accounts 
for the most luminous of quasars and binary X-ray sources. It can also be responsible 
for driving powerful outflows, as we shall sec. Evidence that disks can extend all the 
way down to the central compact object has been provided by ASCA X-ray observations 
of relativistically-broadened fluorescent iron emission lines from low luminosity Seyfert 
and LINER galaxies. (This interpretation has been somewhat controversial; observations, 
with superior sensitivity and spectral resolution, from XMM- Newton are therefore eagerly 
awaited.) 

Not all disks are thin. There are good phcnomenological reasons to suspect that the 
disks contained in many Seyfert and LINER galaxies thicken over some radii to form 
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Figure 1. Broad Fe Ka line from NGC 3516 observed with ASCA by Nandra et aZ.(1999). 



dense, obscuring torii. Theoretically, it has been proposed that the inner regions of disks 
that are supplied with gas at a rate that is either much higher or much lower than the 
Eddington rate (given by MEdd = ^Edd/c 2 = AttGM/ck) will thicken because the gas will 
be unable to cool and the inflow may even become quasi-spherical. Similarly, the disk in 
NGC 4258 is clearly warped and this is thought to be a common occurrence. It has even 
been proposed that radiative torques acting on disks can turn them over (Pringle 1997). 
Polarization observations offer the opportunity to probe the complex geometry of these 
flows. 

Another example of a subtle arrangement of the emitting elements is provided by X- 
ray observations of the thin disks in Seyfert galaxies. The current working model is that 
only a fraction of the binding energy released by the accreting gas is radiated from the 
disk photosphere (with an effective temperature in the keV range). The remainder is 
dissipated in a hot corona, presumably as a consequence of magnetic flaring followed by 
rcconnection and hydromagnetic wave damping. The heated electrons can then scatter 
the escaping soft photons and, as the corona is probably Thomson thick, this will lead to 
a power-law tail in the hard X-ray spectrum. However, roughly half of these photons will 
strike the disk where they can be absorbed if they have low energy and suffer inelastic 
Compton recoil loss at high energy. The reflected spectrum will therefore be convex 
and be imprinted with line features, most famously, the 6.4 keV Ka line of Fe (Fabian 
et a/.2000 and references therein, Fig. (1)). On this basis, it has been argued that the 
widths of lines in Seyfert and LINER galaxies like MCG 6-30-15 imply that the central 
black hole is rapidly spinning. However, the details depend upon the relative location of 
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the region where the Comptonisation is taking place and the region of the disk responsible 
for most of the reflection. There are two ways to test these models using more detailed 
observations. The first is to use a technique called reverberation mapping (e.g. Young 
& Reynolds 2000) . This requires monitoring the variation of the line and the continuum 
simultaneously and deriving the lag in the variation of the former in response to the 
latter. This tells us about the geometry and size of the line-emitting region. The second 
technique is to use polarization observations which will test the geometry 

1.3. Physical Processes 

We now summarize some relevant physical processes. A good general reference is Rybicki 
& Lightman (1979). 

1.3.1. Thomson Scattering 

Classical Thomson scattering is strongly polarizing. A free electron can be considered 
as a Larmor dipole driven by the electric field of the incident wave (with polarization 
vector e). The scattered power into polarization is e* is given by 

— =r 2 (e-e') 2 (11) 

where r e — e 2 /m e c 2 — 2.82 x 10~ 13 cm is the classical electron radius. Note that when 
the scattering angle is <j> = 90°, the radiation is 100 percent polarized. 

Averaging over incident polarization and summing over final polarization states gives 
the familiar differential cross section 

Integrating over solid angle gives the total Thomson cross section 

Q 

ot = —r\ = 6.65 x 10~ 25 cm 2 . (1.3) 

o 

1.3.2. Compton scattering 

At X-ray and 7-ray energy we must take account of the electron recoil. Conserving 
energy and linear momentum, we obtain an expression for the scattered energy e' in 
terms of the incident energy e 

(1.4) 



1 + — ^-r(l -cos< 
Averaging over cj> for small values of e, we obtain the mean energy shift 

e 2 

<Ae>= -. (1.5) 

m e c A 

The Thomson cross section must be replaced by the Klein-Nishina cross section 
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which emphasizes forward over backward scattering. High energy scattering is nearly as 
strongly polarizing as Thomson scattering. Averaging over 4> again for small angle, we 
obtain 

c7 = ct t (i--^2+...) . (1.7) 

V m e c 2 I 
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The rate at which the photons heat the electrons through Compton recoil is therefore 
given for e << m e c 2 by 

/2 
deN(e)-^— = n e o- T cU < £ T , (1.9) 

m e cr m e c^ 

where N(e) is the photon number density per unit energy, U — J deeN is the photon 
energy density and <> should be interpreted as an energy density- weighted photon 
energy. 

These expressions describe the energy loss in the initial rest frame of the scattering 
electron. However when the plasma is hot the electron will be moving and this will cause 
the photon to experience a Doppler shift. To 0(v/c), blue shifts balance redshifts and 
there is no net energy change when the electron distribution is isotropic. However, it is 
apparent what there will be a net energy gain to 0(v/c) 2 = 0(kT/m e c 2 ), because the 
rate of approaching collisions will exceed the rate of receding collisions. We can therefore 
express the net rate of energy loss by the electrons in terms due to the Doppler shift as 

W. = n e o- T cU^—!- , (1.10) 

m e C i 

where a; is a number that we can fix by observing that there should be energy balance, 
W + = W- , when we use a dilute black body of temperature T 

N(e) ex e 2 exp[-e/fcT] . (1.11) 

(It is necessary to use a dilute black body to avoid having to consider nonlinear, induced 
Compton scattering.) As < e >= 4fcT, we deduce that x = 4. 

If this were the dominant physical process, then the equilibrium electron temperature 
in a given radiation field would be 

Tc=^. (1-12) 

However, accretion disk corona are probably heated through reconnection and hydro- 
magnetic turbulence, and the temperature is probably quite non-uniform and hard to 
predict in detail. An additional complication is that the thcrmalization timescales are 
actually quite long compared with the disk dynamical timescale and so the plasma is 
likely to have a significant suprathermal component (Gicrlinski et al. 1999). 

More generally, we can deduce the form of the kinetic equation for the photon dis- 
tribution. As the individual photon energy shifts are small, this will have the form of 
a modified diffusion equation in energy space. However, as the scattering angles are 
typically large, we cannot regard this as a diffusion in momentum space. We there- 
fore just consider an isotropic radiation field to bring out some principles (although this 
approximation is inappropriate for computing polarization, since an isotropic radiation 
field creates zero net polarization upon scattering). As Compton scattering conserves 
the number density of photons the equation must have the form 

dN OF . oN 

7*— &' (L13) 



where F is the flux of photons in energy space. Now for a dilute black body, Eq. (1.11), 
F will be linear in N, as long as we can ignore induced scattering, and, as it represents 
a diffusion, F will contain the first derivative of N with respect to energy. As F must 
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vanish for a dilute black body it must have the form 

for some function g{e). We next multiply this equation by e and integrate over energy 
and use either Eq. (|l.9[) or Eq. (1.10) to identify the function g(e) = n e aTce 4 kT/m e c 2 . 



The resulting (Fokker-Planck) equation is, 

/fXNl^W 1 

(1.15) 



dN n e UT d 2 
dt m e c de 



e»*r(^l + .v 



This formalism describes the behavior of electrons interacting with a dilute gas of photons 
inside a box with reflecting walls. If we need to take account of induced scattering 



processes, then Eq. ( 1.14 ) must be modified soF = when N(e) has the Planck form, 

AT( e ) = |^![exp(e/feT)-l]- 1 . (1.16) 



The result is that N — > N(l + N), in the second term in brackets in Eq. (1.15) which is 
then known as the Kompaneets equation and is central to discussions of the transfer of 
radiation through hot plasmas. 

Clearly the radiation as described by this isotropic formalism will be unpolarized. 
In order to describe the polarization of a cosmic source, we must tackle the radiative 
transfer. There are three approaches that have commonly been followed. 

(a) Escape Probability F orma lism This is the simplest approach. We add a term to 
the right hand side of Eq. ( 1.15J ) 

Nc 

-wh- ,L17) 

The extra factor 1 + t takes into account the impeded photon escape when the Thomson 
depth t exceeds unity. This approach, which is most commonly used, although instructive 
as far as the spectrum goes, is not much help when it comes to polarization. 

(&) Intensity Formalism Provided that we restrict our attention to simple shapes - 
slabs, spheres etc. - we can incorporate the energy space transport within the equation 
of radiative transfer using a scattering kernel. This can then be solved by taking mo- 
ments and imposing a closure relation in the standard manner. It is possible to include 
polarization though this leads to quite involved equations. 

(c) Monte Carlo Formalism In many ways the most versatile method is the same 
one used in nuclear reactors for the transport of neutrons. This is to follow individual 
photons, within the scattering region starting with energies, locations and directions that 
are selected according to a prescribed distribution using random numbers. Polarization 
is relatively easy to handle, and most polarization is created in the last few scatterings, 
reducing the computational burden which occurs at large optical depths. All of this is 
quite straightforward in principle, though, in practice, Monte Carlo simulations are quite 
an art as a variety of ingenious tricks have to be used to reduce the variance with a finite 
amount of computer time. 

1.3.3. Dust Scattering 

Dust scattering is more complex than Thomson scattering as the cross section depends 
on wavelength, grain size, and grain composition (cf Hildebrand, these proceedings). In 
the limit when the wavelength is much greater than the grain size, Rayleigh scattering 
applies, which has the same angular cross section as Thomson scattering, but scales as 
A -4 , a fact which has been used to distinguish electron scattering from dust scattering in 
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some Seyfert 2 galaxies. For a range of different sizes and compositions, it is impossible 
to express the dust scattering cross section in a simple formula, but extensive numerical 
calculations have been carried out by, e.g., Draine & Lee (1984), Zubko & Laor (2000). 

1.3.4. Faraday Rotation 

The next relevant physical process is Faraday rotation. When electromagnetic wave 
modes propagate through a plasma, their phase velocities will be changed from c. To 
lowest order, the eigenmodes are circularly polarized with phase velocity difference 

AT4 = 2c-^cosa, (1.18) 

where lo p — (47rn e e 2 /m e ) 1 ' 2 is the electron plasma frequency, ujq is the electron gyro 
frequency and a is the angle between the ray and the magnetostatic field. If we decompose 
a linear polarized wave into two circularly polarized modes that propagate through the 
medium and then recombine the modes after they leave the medium, then there will be 
a net rotation of the plane of polarization through an angle 

Polarization observations can tell us as much about the intervening medium as about 
the source. 

1.3.5. Relativistic Radiative Transfer 

Another interesting complication is that special and general relativistic effects will 
affect the transfer of radiation from the disk to us. This is particularly interesting for 
line radiation. The non-relativistic Doppler shift will broaden the profile of a line formed 
at the photosphere of a rotating disk, with the blue wing coming from the approaching 
limb and the red wing from the receding limb. The gravitational shift (which is not 
separated from the Doppler shift in a general relativistic calculation) will accentuate the 
red wing whose extent depends upon how close the inner edge of the disk gets to the 
hole. A further effect is that rays will be deflected by the gravitational pull of the central 
black hole so that an image of an accretion disk observed from near the equatorial plane 
would exhibit the back side of the disk. (There are ambitious proposals to deploy an 
X-ray interferometer in space which could exhibit these and other effects.) 

Of more direct relevance to this school, is the behavior of the polarization. There 
are two main effects. Firstly, aberration changes the emission angle, and consequently 
the emitted polarization, in the rest frame of the orbiting gas. Secondly, the plane of 
polarization is rotated as the ray propagates near the black hole (Fig. (2)). As we discuss 
in more detail in §5, the polarization direction is parallel-propagated along null geodesies. 
All of this is straightforward, if somewhat tedious, to compute. Specific models for 
continuum emission have been computed by Laor et al. (1990) and Connors et al. (1980), 
while relativistic effects on line polarization were computed by Chen & Eardley (1991). 
It may be possible to use observations of the rotation of the polarization direction with 
wavelength, in a spectral line, or in the continuum if the wavelength is a measure of the 
effective radius of the disk, to measure the spin of the hole. 

1.4. Interpretation 

Having outlined some of the relevant physical mechanisms, let us return to the problem 
with which we began this section. "How much polarization do we expect from an accre- 
tion disk and what do we observe?" We can split the problem into two parts - the escape 
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Figure 2. Polarization angle as a function of position for an electron-scattering dominated thin 
accretion disk around a Kerr black hole (a=0.998) viewed at infinity from an angle of 75°. The 
figure is 20GM/c 2 in size. 



of photons emitted in the disk and the behavior of the scattered photons. For a given 
atmosphere, the problem is linear and we can superpose the two components. However, 
if we try to solve self-consistcntly for the ionization and thermal state of the atmosphere, 
the problem becomes nonlinear. 
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1.4.1. Electron Scattering 

Radiative transfer in a pure scattering, plane parallel atmosphere is a classical problem 
that was solved for a Thomson scattering kernel, initially analytically, by Chandrasckhar 
(1960) and then in greater generality by Angel (1969) using a Monte Carlo approach. 
The answer is that the emergent polarization varies with inclination, having a value 
p = 0.12 when the atmosphere is viewed horizontally and p — 0.02 when viewed at the 
most probable inclination of 60° and, of course, p — when viewed normally. Real 
disk atmospheres also have an absorptive opacity and this will reduce (or increase) the 
emergent polarization significantly (Hubcny et aZ.2000). In addition, we now believe that 
accretion disks are strongly magnetized. The rationale for this is that ionized accretion 
disks are known to be unstable to developing strong internal magnetic fields with interior 
magnetic pressures estimated to be ~ 1 — 10 percent of the gas pressure. This magnetic 
field will surely be carried out beyond the photosphere and into the coronae (discussed 
above) by buoyancy forces. Furthermore, magnetic pressure is likely to dominate gas 
pressure in an accretion disk corona, just like in the solar corona. 

When we consider the specific parameters appropriate to observed disks, we find 
that thermal emission should be unpolarized based on the following argument. Con- 
sider a given disk annulus, the radiation pressure at the photosphere, aT 4 /3, is smaller 
than the gas pressure within the disk. Because observations at a given thermal wave- 
length peak near A ~ hc/kT, we can derive a lower limit on the magnetic field strength 
B ;> 10 2 (A/500nm)~ 2 G. Since photons traverse differing optical depths after their last 
scattering and the magnetic field will likely have significant inhomogencities, any po- 
larization will be erased by the Faraday rotation (d<f> /oItt) J> 10, independent of the 
wavelength observed or the physical size of the disk (Agol & Blaes 1996). 

Thus, it should come as no surprise that AGN disks are generally only polarized by a 
small amount (Antonucci, this volume). Furthermore, the polarization that is observed 
may be imprinted extrinsically. Purported rises in polarization below 912 A in a few 
quasars observed with HST contradict old theoretical predictions and stand as a challenge 
to disk theory (Koratkar & Blaes 1999). 

Turning to X-ray wavelengths, where non-thermal emission means that we might ignore 
Faraday rotation, but cannot ignore reflection, predictions of the emergent polarization 
under a variety of models have been presented by Matt, Fabian, & Ross (1993) and 
Poutanen & Svensson (1996) and references therein, as shown in Fig. (3). The rather flat 
X-ray spectrum is created by Compton scattering of thermal emission from the accretion 
disk, which is partly due to absorbed X-rays. The electron-scattering reflection feature 
is suppressed at low energies, <; 8 kcV, by X-ray bound-free absorption opacity, and 
at high energies, ^100 keV, reduced by electron recoil; consequently, the largest X- 
ray polarizations should lie between these energies. The magnitude of the polarization 
may depend upon the exact geometry and placement of the coronal emission regions, an 
unexplored problem. 

1.4.2. Dusty Disks 

The outer parts of accretion disks may be cool enough (T < 1800 K) to be inhabited 
by dust grains. If the dusty disk drives a wind, or is inflated or warped, then dust will 
scatter the light from the inner disk, imprinting a polarization signature from the infrared 
to ultraviolet. It may in practice be quite difficult to distinguish between a dusty disk, 
torus, or outflow using polarization, as the level of polarization expected is quite small, 
~ 1%, and depends on the details of the dust model, e.g. Konigl & Kartje (1994). In 
addition, dust extinction can create polarization if the grains, charged by collisions with 
ions, arc aligned by magnetic fields, inducing polarization at the percent level as well. 
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Figure 3. Spectrum and polarization versus frequency (x = hv/m e c 2 ) of hot, plane-parallel 
corona (T e = 0.11 m e c 2 , tt — 0.05) above a cold disk viewed at two inclination angles, 
H — cosi = 0.11 and 0.5. The numbers label various scattering orders for reflection, while 
the solid lines show the total spectrum and polarization (from Poutanen & Svensson 1996). No 
relativistic propagation effects have been included. 

1.5. Summary 

• Disks are commonly found in accreting systems. 

• Model accretion disks can create strong polarization both in transmission and in 
reflection, throughout the electromagnetic spectrum. 

• Measurement of the variation of linear polarization with wavelength can, in principle, 
reveal a lot about the disk structure and the location of coronal emission sites. 

• However, the situation is, in practice, more complex, particularly at optical wave- 
lengths, where external illumination, warping, and especially Faraday rotation are likely 
to be very important. 

• Even when polarization cannot be measured, its effects are so large under conditions 
of strong electron scattering that radiative transfer calculations should include polariza- 
tion. 

• Monte Carlo techniques are well-suited for computing polarization in a given model. 

• There is a very strong case for developing X-ray polarimctry. 



2. Jets 



2.1. Motivation 



Jets, or more generally bipolar outflows, are also surprisingly common. They have been 
studied in association with active galactic nuclei (AGN), binary X-ray sources, young 
stellar objects, novae and so on. Jet speeds are typically a few times the escape velocity 



Roger Blandford et al.\ Compact Objects and Accretion Disks 11 

from the central object; in the case of black holes, bulk Lorentz factors of 7 ~ 10 are 
inferred. (Gamma ray bursts may also produce jets with Lorentz factors 7 ~ 300.) Jets 
are so common that it has been speculated that they may be an essential concomitant 
of accretion flow - the channels through which the liberated angular momentum and 
perhaps also much of the energy leave the system. The challenge to the astrophysicist 
is to explain how jets are powered and collimated. However, even after decades of work, 
major theoretical and observational questions about their origin, collimation and even 
their constituents still remain. 

There are, generically, two proposed origins for the jet power: the central object (black 
hole, neutron star or protostar) and the accretion disk (e.g. Blandford et al.1990). In 
both cases, the energy derives from differential rotation. For example in the case of a 
Keplerian disk that extends down to the surface of a non-rotating, unmagnctized star, as 
much energy is released in the boundary layer as in the disk. An extreme case is presented 
by the Crab (Weisskopf et aZ.2000) and Vela (Helfand, Gotthclf & Halpern 2001) pulsars 
which exhibit prominent jets without there being any accretion disk, presumably. 

One of the best observational approaches to investigate the mechanisms which produce 
jets is to determine the jet composition at radii where they can be observed directly. In 
the case of black-hole jets, the plasma is likely to be electron-ion if the jet originates 
from a disk, or electron-positron if it derives from the black hole. Polarization obser- 
vations have been prominent in attempts to distinguish between these two possibilities. 
Most contemporary explanations of the collimation invoke strong magnetic fields, though 
in most cases, the argument for magnetic collimation is a theoretical one, based upon 
eliminating the alternatives. One exception to this is the bipolar outflow associated with 
young stellar objects (YSOs) where polarization observations strongly support the notion 
that the magnetic field is dynamically important (Akeson & Carlstrom 1999). 

2.2. Observations 

The first jet observed was that in the nucleus of the elliptical galaxy M87 in the Virgo 
cluster. The modern representation is given in Fig. (4) (Perlman et aL1999). In this 
case, the jet emerges from no more than 100 times the gravitational radius of the central 
hole (M = 3 x 10 9 M Q ; GM/c 2 = 4 x 10 14 cm) and propagates radially outward for a 
distance ~3x 10 22 cm, seven orders of magnitude larger. The jet, however, is neither 
homogeneous nor smooth. Superimposed upon an overall decrease in surface brightness 
as the jets expand away from the central hole are strong side to side variations and 
bright features. These may reflect a time-dependence at the jet origin or independent, 
local instabilities. At small radii, the M87 jet is strongly one-sided, and this is generally 
attributed to Doppler beaming and adduced as evidence for relativistic outflow. 

Another good example of an extragalactic jet is Pictor A, which has recently been 
observed at X-ray energies using Chandra (Fig. (5)); Wilson, Young & Shopbell, 2001). 
It shows strikingly efficient collimation and powerful emission from the western hot spot, 
which is also prominent at optical and radio wavelengths. The Galactic source GRS 
1915+105 shows mild supcrluminal motion which is interpreted in terms of moving fea- 
tures, probably internal shocks, with space velocities ~ 0.9c (e.g. Mirabel & Rodriguez 
1994). This is mild compared with most extragalactic compact jets where the speeds 
are more typically ~ 0.99 c, corresponding to bulk Lorentz factors 7 ~ 7. Indeed, there 
are observational indications that much larger bulk Lorentz factors are produced in ex- 
tragalactic jets. Independent evidence that these jets are relativistic comes from 7-ray 
observations, principally done with the EGRET detector on Compton Gamma Ray Ob- 
servatory. These showed that those jets that are beamed towards us, and which are 
collectively known as blazars, are often powerful 7-ray sources. The total electromag- 
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Figure 4. False-color representations of the total intensity and polarization of the M87 Jet 
in the optical (HST F555W, top two panels) and radio (VLA 14.5 GHz, bottom two panels). 
The HST observations were carried out in May 1995, while the VLA observations were done in 
February 1994. All maps were rotated so that the jet is along the x-axis, and are convolved to 
0.23" resolution (from Perlman et al. 1999). 



netic spectrum of blazars (and similar sources) comprises a broad band synchrotron 
radiation spectrum extending from low radio frequencies to an upper frequency between 
optical and X-ray wavelengths. The same electrons are responsible for an inverse Comp- 
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Figure 5. Gray-scale representation of the full-resolution Chandra image of the nucleus, jet 
and western hot-spot of Pictor A (from Wilson et a/.2001). 

ton component that can extend up to TeV energies (with variability times as short as 
~ 30 min). (We can only observe TeV emission from relatively local sources because TeV 
photons from cosmologically distant sources will be absorbed through pair production on 
the intergalactic infrared background.) The fastest jets may well be associated with 7-ray 
bursts, if they are indeed beamed, for which speeds of ~ 0.999995c have been inferred. 

2.3. Physical Processes 

2.3.1. Synchrotron Radiation 

As described in Dr. Landi Dcgli'Innoccnti's contribution (for more details see also Ry- 
bicki & Lightman 1979), the synchrotron power radiated by an ultra-relativistic electron 
with energy jm e c 2 in a field of strength B is given by 

P = ^ 2 <7 T cUb , (2.20) 

where Ub = B 2 /8ir is the magnetic energy density and we have averaged over pitch 
angle. The corresponding radiative cooling time is 

U-'A, (2.21) 

where B is measured in Gauss. The characteristic frequency radiated is 

i/ c = 7 2 BMHz. (2.22) 

The polarization of single particle emission varies from 2/3 for v <C v c to 1 for v 3> v c . 
The electric vector is perpendicular to the projected magnetic field direction. Averaging 
over a power law distribution of relativistic electrons, dN/d'y = K j~ s , it can be easily 
shown that the observed intensity is 

I„xKB 1+a i>- a , (2.23) 

where the spectral index a = (s — l)/2. The net degree of linear polarization is 

Q 8 + 1 



I s + 7/3 



(2.24) 
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Synchrotron radiation is also naturally circular polarized to an extent 

j ~3/ 7) (2.25) 

dependent upon the detailed angular distribution function and the viewing angle. 

An important issue for what follows is the viability of a synchrotron maser. This can 
be shown to be impossible for ultra-relativistic emission in vacuo. Essentially, in order 
for a maser to operate, it is necessary that there be a population inversion and that the 
emissivity at a given frequency decreases sufficiently rapidly with increasing energy and 
this does not happen with regular synchrotron emission. (It can however arise when 
relativistic electrons emit synchrotron radiation in a plasma, although in practice the 
conditions for this to occur are rather restrictive. This is known as the Razin effect.) 

Of more relevance is what happens when the electrons are no more than mildly rela- 
tivistic. The emission is then confined to a series of harmonics of the fundamental gyro 
frequency u>g/j- As the central frequency of a harmonic decreases with increasing energy, 
there are frequencies and directions where maser action is possible, given a population 
inversion. Cyclotron masers are likely to be highly polarized. The polarization from 
an electron-ion plasma will be elliptical; that from an equal pair plasma will be purely 
linear. 

When the brightness temperature Tb = I v c 2 /2kv 2 of a synchrotron source approaches 
the kinetic temperature of the emitting electrons, T^ — r ym e c 2 /3k oc (ly/B) 1 / 2 , the 
radiation will be absorbed. The optically thick radiation from a source will have a 
brightness temperature that is limited to this value and so the optically thick equivalent 
of Eq. (HH) is 

J„ oc J o/2 B- 1/2 . (2.26) 

When we consider the linear polarization of a self-absorbed source, we observe that the 
electrons emitting in the field-perpendicular polarization at a given frequency will have 
slightly lower energies than those emitting in the field-parallel polarization. Therefore 
the brightness temperature of the field parallel emission will be slightly larger than that 
of the field perpendicular emission. The degree of linear polarization from a power-law 
relativistic electron distribution function can be computed to be 

Q -3 



I 6s + 13 



(2.27) 



An example is the supermassive black hole candidate in the Galactic center, Sagittar- 
ius A* , which has no linear polarization up to frequencies of 86 GHz (cf Hildebrand, these 
proceedings), but shows surprisingly strong circular polarization (e.g. Bower 2000). The 
degree of CP increases sharply with frequency. Whereas both advection-dominated- 
accretion-flow (ADAF) models and accretion-disk-powered-jet models can account for 
the spectrum from centimeter wavelengths to X-rays of Sagittarius A*, they have dis- 
tinct polarization characteristics which might in the near future be able to distinguish 
between the two models. 

2.3.2. Inverse Compton Scattering 

Inverse Compton (or more properly Thomson) scattering, in which a highly energetic 
particle transfers momentum to a low-energy photon, is very similar to synchrotron 



radiation. If the radiation field is isotropic, then the power is given by Eq. (2.20) with 
the magnetic energy density Ub replaced by the radiation energy density, £/ ra d. As 
photons are conserved in Thomson scattering, the mean photon frequency is boosted by 
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an average factor 

- = \l\ (2.28) 

v 3 

where v' is the scattered frequency and v is the incident frequency. One power of 7 arises 

from the Lorentz transformation into the electron rest frame; the second comes from the 

scattering back into the original frame. The polarization observed will be generically 

be ~ I/7 unless the incident radiation field is both highly anisotropic and polarized, in 

which case a strongly polarized scattered spectrum can be emitted. This can arise when, 

for example, radiation is scattered into a beam by a warped disk. Similarly, circular 

polarization in the incident radiation will be partly retained in the scattered radiation. 

2.3.3. Inverse Compton Limit 

The comparison of synchrotron radiation and inverse Compton scattering leads to 
what, for historical reasons, is called the inverse Compton limit (Kcllermann & Pauliny- 
Toth 1969). The way the argument is traditionally expressed is that the ratio of the 
Compton power radiated by an electron to the synchrotron power can be written as 



L c -i _ Us_ Vs T ' 
L s ~ Ub^ B 2 



<x -±-f- oc T b v s , (2.29) 



in obvious notation, and where we have assumed that the source is self-absorbed at the 
observing frequency. If wc set v$ ~ 1 — 10 GHz, then the brightness temperature is 
limited to Tg ~2x 10 12 K if this ratio is not to exceed unity. 

The original concern was that if the ratio did exceed unity then the second order 
Compton scattering would be even greater than the first order scattering and so on. Of 
course this can't go on for too many orders because the Klein-Nishina limit will limit the 
scattering. Nowadays we think we can identify the synchrotron and the inverse Compton 
components and so we know their ratio and can deduce the source brightness temperature 
which is quite insensitive to its value. 

2.3.4. Kinematics of Bulk Relativistic Motion 

In the case of a relativistic jet, it is often easier to compute the radiation spectrum in 
the comoving (primed) frame of the emitting plasma and then perform a Lorentz boost 
into the (unprimed) frame of the observer. The frequency will be boosted by the Dopplcr 
factor (e.g. Blandford & Konigl 1979), 

where 9 is the scattering angle in the observer frame and (3 = v/c is the bulk velocity of 
the plasma in the jet. Note that for a jet beamed toward us with 9 < 7 _1 , the Dopplcr 
factor is <5 ~ 7. Note also that the rate of change of observer time t bs to proper time r 
satisfies 

£- = 5. (2.31) 

As a consequence, the observed transverse speed of a feature moving with the jet speed 
is given by 

/U s = 7/3sin^ = -^fli-. (2.32) 

creobs 1 — p cos 9 

This has a maximum value 7 (3 for 9 — cos _1 /3 and consequently the expansion can 
be "supcrluminal" when (3 > 0.71. The kinematics of real jets is undoubtedly more 
complex and the space motion of shock features must be distinguished from the speed of 
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Figure 6. The Lorentz boost of polarization. 



the emitting plasma. Frequently, observers make the approximate identification /3 bs ~ 
9~ x ~ 5 ~ (dt/dtobs) ~ 7 when interpreting measurements of compact extragalactic 
radio sources. 

The behavior of polarization under a Lorentz boost is straightforward. The k-vector 
swings forward along the direction of motion making an angle 9' with the direction of 
motion in the plasma rest frame and an angle 9 in the jet frame, where sin 9' = S sin 9. If 
we imagine the k-vector as being rotated in this manner, then the electric and magnetic 
fields associated with individual photons will be similarly rotated about a direction kx B 
so that k, E and B continue to form an orthogonal triad (Fig. (6)). 

Because the Planck distribution function is Lorentz scalar, and the brightness tem- 
perature enters only in the ratio v /T, it is clear that the brightness temperature must 
transform in the same manner as frequency. This implies that if the inverse Compton 
limit is applied in the observer frame, the brightness temperature measured by a radio 
astronomer can be as high as ~ 2 x 10 12 5 K. This is particularly germane at this time 
because radio astronomers are able to estimate these brightness temperatures, both di- 
rectly using ground and orbiting VLBI, and indirectly by carefully analyzing refractive 
interstellar scintillation. Using the limits on the source size derived from these obser- 
vations and its relation to the source flux-density (i.e. the Rayleigh- Jeans equation), 
one can derive an apparent surface brightness temperature and therefore the value of 5, 
assuming the intrinsic brightness temperature cannot exceed the inverse Compton limit 
significantly. The values that are found require bulk Lorentz factors 7 ~ 5 ~ 30 and may 
be even higher. 
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Figure 7. A superposition of self-absorbed sources peaking at different frequencies can create 

a combined 'flat' spectrum. 

2.3.5. Faraday Conversion 

We discussed Faraday rotation in a cold, non-relativistic plasma in the last section. 
We must now consider what happens in an ultra-relativistic pair plasma. On symmetry 
grounds, the eigenmodes must be linearly polarized and are usually labeled ordinary, 
where the electric vector along the direction k x B, and extraordinary, where it is not. 
If a linearly polarized wave, obliquely polarized with respect to the magnetic field is 
incident upon the plasma, then it can be decomposed into ordinary and extraordinary 
modes that will propagate with slightly different phase velocities. In this manner, circular 
polarization will be created. In other words, there is a conversion of U to V. The sense 
of circular polarization will be given by the sign of (e • B){e ■ k x B) and so in order to 
have a measurable circular polarization from a cosmic source, it is necessary to have a 
preferred field orientation. This phenomenon is known as Faraday conversion (e.g. Jones 
& O'Dell 1977). 

Because Faraday conversion is caused by the lowest energy relativistic electrons, it 
can serve as a probe of the low energy end of the electron energy distribution. Faraday 
conversion is furthermore proportional to e 2 B 2 /m 2 , i.e. independent from the sign of 
the particle's charge. An equal mixture of electrons and positrons can therefore produce 
Faraday conversion, but not rotation. A comparison of linear and circular polarization 
might therefore probe the constituents of the jet (e.g. electron-positron versus electron- 
ion pairs). 

2.4. Interpretation 

2.4.1. Shocks and the Integrated Spectrum 

The emitting element in powerful synchrotron jets is thought to be a relativistic shock 
wave and a typical source will comprise the emission from several of them. In the limit, 
a jet can be thought of as accelerating relativistic electrons over its length with a field 
strength that diminishes with radius. The total flux density observed at a given frequency 
is dominated by the emission from the radio photosphere, where the optical depth is 
roughly unity. Blazars and similar sources generally have "flat" spectra, i.e. —0.5 < a < 
0.5 and can be interpreted as the superposition of a series of self-absorbed sources at 
successive radii, each peaking at successively lower frequencies (Fig. (7)). 
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2.4.2. Jet Composition 

On the basis of the fraction of the jet energy per relativistic electron responsible for 
the synchrotron radio emission that we observe to be radiated on an outflow timescale, 
it has been tentatively deduced that jets cannot carry protonic "baggage" and so must 
comprise electron-positron pairs with a low energy cutoff in the distribution function 
(e.g. Reynolds et aL1996). 

Wardle et aZ.(1998) measure a large degree of circular polarization in 3C279, a source 
where the linear polarization is also quite high, and so there cannot be too much normal 
Faraday rotation. In order to explain the circular polarization, they have to invoke a 
large population of mildly relativistic electrons and positrons. On this basis they conclude 
that relativistic jets comprise pair plasmas at least at the radii where they are directly 
observed. 

2.4.3. Coherent Emission Mechanisms 

All of this calls into question the fundamental synchrotron hypothesis. If the deduced 
jet powers are unreasonably large then we should certainly be prepared to consider the 
possibility that the radio emission, or at least its compact and variable part, may be 
due to a coherent emission mechanism. This is not unreasonable. After all, the sun 
and Jupiter support high brightness coherent emission under far more docile conditions. 
Furthermore, a shock front is a very natural environment in which strong and unstable 
currents are likely to be induced and these are commonly observed to radiate coherently 
in plasmas. Probably the most likely possibility is coherent cyclotron emission emerging 
from very much more compact regions than under the jet hypothesis. This requires the 
magnetic field strength to be hundreds or even thousands of Gauss and so the pressures 
must be much larger than in the jets. Coherent cyclotron emission is likely to be strongly 
circularly polarized, unless there are equal numbers of electrons and positrons. 

2.4.4. Microlensing and Refractive Scintillation 

Both microlensing by compact objects in the line-of-sight and refractive scintillation by 
density fluctuations in the Galactic ionized medium can introduce non-intrinsic variability 
of the compact structures in jets of extra-galactic radio sources (e.g. knots or shock 
fronts). For significant variability to occur, these structures must have an angular size of 
the order the Fresnel scale for scintillation or Einstein radius in the case of microlensing 
(e.g. Koopmans & de Bruyn 2000). Both scales are typically several micro-arcseconds. 
The expected time-scale of variability is determined by the transverse velocity of the 
source compared to the scattering medium (i.e. the compact objects or the Galactic 
ionized ISM) and is typically hours to weeks in the case of refractive scintillation or 
weeks to months in the case of microlensing. 

Whereas polarization is typically little affected by either scintillation or microlensing 
(which retains polarization angle), in both cases only the most compact source struc- 
tures vary significantly. If these structures have different polarization degrees or angles 
compared with the flux-density weighted average over the source, the net results will 
be polarization variability that strongly correlates with the non-intrinsic source varia- 
tions. A correlation between changes in polarization and non-intrinsic flux variations of 
extra-galactic radio sources (for example in intra-day variables (IDVs), which strongly 
scintillate) could therefore provide information on the polarization properties of the most 
compact micro-arcsecond scale jet structures, which are impossible to observe directly in 
any other way. 

There are several other ways that the Galactic ionized ISM can introduce non-intrinsic 
variations in the polarization of radio sources (including radio jets), or even induce cir- 
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cular polarization. The simplest case is that of extreme scattering events (ESEs), where 
large localized overdensities in the Galactic ionized medium move into the line-of-sight to 
the radio source. The enhanced electron column density increases the Faraday rotation 
and could result in an observable change in the polarization angle of the radio source, as 
well as in a change of its flux-density due to refractive lensing. Besides Faraday conver- 
sion, which converts linear polarization to circular polarization, if a strong gradient in the 
rotation measure over the source exist, the scintillation patterns of left and right-hand po- 
larized wavefronts will be slightly displaced. This will introduce a time-variable circular 
polarization that is strongly correlated with the scintillation-induced flux-density varia- 
tions (Macquart & Melrose 2000), but independent of the degree of linear polarization 
(which is not the case for Faraday conversion). 

2.4.5. Inverse Compton Scattering 

A major concern in interpreting the inverse Compton X-ray and 7-ray observations 
of blazars is the source of the incident photons. In the lower power objects that can be 
observed at TeV energies, these are thought to be synchrotron photons emitted locally 
within the jet. In the higher power quasars where there is a powerful photoionizing 
continuum as well as a relativistic jet coming towards us, most of the incident photons 
are thought to originate from the disk and to be scattered into the jet. However, this is 
not certain and X-ray (or 7-ray) polarimetry could validate this, because these scattered 
X-rays of external origin would be highly polarized. 

2.5. Summary 

• Jets appear to be a common and perhaps even a necessary feature of accreting or 
possibly even just simply rotating systems. 

• Despite much observational progress, the fundamental questions concerning the ori- 
gin, composition and collimation of jets remain unanswered. 

• Black hole jets are formed ultra-rclativistically as required to account for their su- 
perluminal motion, high radio brightness temperature and prodigious 7-ray emission. 

• Radio polarimetry of relativistic jets is helping us to deduce their composition al- 
though coherent emission mechanisms cannot be ruled out. 

• X-ray and 7-ray polarimetry could help refine models of relativistic jets by probing 
jets close to their origins. 



3. Outflows 

3.1. Motivation 
Somewhat paradoxically, accreting systems commonly exhibit outflows. The reason for 
this behavior is simple. As gas accretes onto a compact object it must release its grav- 
itational binding energy. When this is possible, it will do so by radiating. However, 
this may not be possible when the gas accretes much faster than the Eddington rate, 
A^Edd = 4:1tGM/ktc, the photons will be trapped by electron scattering and the energy 
can only be carried off by a bulk outflow. Even if the radiation is not trapped, then 
it can still drive an outflow if it encounters gas with an opacity k >> kt, for example 
with dust grains or resonance lines. Either, or more likely both, of these processes are 
believed to drive the outflows associated with Broad Absorption Line Quasars (BALQs) 
and this lecture will be primarily about these objects, although the principles involved 
are more generally applicable. The broad absorption lines by which these quasars are 
distinguished are associated with the ultraviolet resonance lines of the common ions. 
(Recall that quasars are mostly at high redshift and so these lines are observed in the 



20 



Roger Blandford et al.: Compact Objects and Accretion Disks 



• NELR 




Seyfert 2 



Figure 8. Orientation model for Type 1 and Type 2 Sefert galaxies. In the simplest scheme, 
only the narrow line region can be seen from a Setfert 2 galaxy; the broad line region lies behind 
a thick equatorial ring of obscuring gas and dust. 



visible.) They show absorption troughs, extending to the blue of the regular, broad 
emission lines by which quasars arc identified spectroscopically. This is just what is also 
seen in the star P Cygni, though the BALQ relative velocities, ~ 0.1c, which are what 
one might escape for gas escaping from the vicinity of a black hole, are much larger than 
those encountered in P Cygni (Weymann et al.1991). 

However, we do not understand the flow of gas around the black hole and the location 
of the emission line and the absorption line gas is still quite controversial. Models of broad 
line clouds have been constructed in which the gas flows "in, out, round or about" and 
there may be elements of each of these four kinematic classes in real sources. It is probably 
safe to conclude that the flow is not very simple; otherwise we would have already 
understood it using a technique called "reverberation mapping". In this technique it 
is supposed that the lines vary in direct response to the photoionizing continuum. By 
monitoring them both, it is possible to construct a Greens function response of the 
emission line gas, This can then be compared with the predictions of simple kinematic 
models. It has been possible to use this technique to locate the gas in several instances. 
However, the details of the velocity field remain controversial. 

Despite this, most authors have assumed that the outflow of the absorbing gas is 
roughly equatorial and radial, originating from close to the hole. This implies that the 
quasar will only be classified as a BALQ when the observer direction is also close to the 
equatorial plane; otherwise she will see a regular radio-quiet quasar (Fig. (8)). The goal 
is to use polarization observations to see if this is truly the case (Antonucci 1993). 
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3.2. Observation 

As we have remarked, a typical BALQ spectrum shows broad emission lines like CIV 
A1548, accompanied by broad troughs extending to the short wavelength end of the 
spectrum with widths up to ~ 30,000 km s _1 . The emission lines are generally unpo- 
larized, though both the semi- forbidden line CIIIJA1909 and Lya can be polarized, as 
we shall discuss. The troughs which arc caused by approaching gas, represent photons 
that are removed from the radiation field and scattered sideways. They also remind us 
that momentum is taken out of the radiation field so that the gas is accelerated (Arav 
et aL1995). (Of course this may not be the only accelerating force involved, though it is 
simplest to assume that it is.) 

The absorption troughs themselves are not black which presumably means that gas 
moving along different directions is scattering radiation into our line of sight. In addition, 
although the troughs are relatively smooth, they do show velocity structure and a variable 
degree of polarization that can be as large as ~ 0.2 a characteristic feature of resonance 
scattering. However, resonance scattering is not the only means of producing linear 
polarization; electron scattering can do the job just as well and some observers have 
preferred this explanation. 

In what follows, we shall confine our attention to resonance scattering for two reasons. 
Firstly, we observe resonance scattering directly and its contribution to the opacity is 
three to four orders of magnitude larger than that of electron scattering. Secondly, 
the physics is much more interesting than that of electron scattering! Furthermore we 
shall confine our attention to the polarization and the implications that it has for the 
kinematics, as opposed to the dynamics of the flow. 

However, we must mention some additional observational clues as to the nature of 
BALQs. Firstly, BALQs are both radio- and X-ray-quict. In some models this implies 
that there is a highly ionized region that can absorb the X-rays and transmit the ul- 
traviolet radiation. In other interpretations, the X-rays are never emitted in the first 
place. Secondly the BAL phenomenon is pretty much confined to quasars; the lower 
power Seyfert galaxies do not exhibit these broad troughs. Thirdly, the continuum (after 
subtracting the galactic contribution) is fairly uniformly polarized with p ~ 0.01 — 0.05, 
suggesting electron scattering in a disk corona is responsible. This is presumably located 
inside the absorption line gas. As mentioned above, it is also possible that this same 
scattering occurs outside the absorption line gas in which case it would probably be re- 
sponsible for the polarization of the absorption line troughs. Finally, a large polarization, 
increasing towards shorter wavelength, has been reported at wavelengths longward of the 
Ly edge in a few, regular quasars. (It should be emphasized, though, that these HST 
observations were very difficult to make and the results have been controversial.) 

3.3. Physical Processes 
3.3.1. Resonance Transitions 

Let us first consider the levels of intermediate Z ions adopting the Russell- Saunders 
approximation. We can distinguish the Li-like ions, such as CIV, NV, OVI, the Be-like 
ions, such as CIII, NIV, OV and the B-like ions such as CII, NIII, OIV. The electronic 
state is determined by the principal quantum numbers for the valence electrons (n,£), 
the total orbital and spin angular momenta (L, S) for the term and the total angular 
momentum J for the level which, in turn is divided into 2 J+l sublevels. This information 
is encoded in the quantum mechanical designation of the level. For example, the ground 
state of OIV is 2s 2 2p 2 P°/ 2 which means that one of the three valence electrons is 
unpaired in a 2p level with total orbital angular momentum L — 1 (i.e. a P state) and 
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DOUBLET J=l/2 to 1/2, 3/2; eg CIV 1558, Ly 0C1216 
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Figure 9. Atomic transitions that contribute to the overall polarization for a typical dou- 
ble transition. If the doublet is resolved then the individual components have perpendicular 
polarizations of 3/7 and 0. If it is not, then the average polarization is 3/11 



total spin S = 1/2, hence the superscript 25 + 1 = 2. There are two possible choices 
for J and the lower energy one, by Hund's rule, has J = 1/2 as designated in the final 
subscript. (The final superscript, o, indicates an odd parity, which must change under 
a permitted (electric dipole) transition.) The electric dipole selection rules are that 
A£ = ±1; AS = 0; AL = 0, ±1; A J = 0, ±1, (except that an J = ->• transition is 
forbidden); and AM — 0, ±1, (except that a M = — > transition is forbidden). 

3.3.2. Singlets 

Now consider a singlet transition such as the Be-like CIIIA977. The ground state 
has J = and only one sublevel with M = 0; the excited state has J = 1 and only 
permitted transitions to M — ±1 need be considered. The radiation pattern depends 
only upon the angular parts of that wavefunction, through the Wigner-Eckart theorem. 
When the scattering angle 9 = n/2, the degree of linear polarization can be computed 
to be p(n/2) — 1. For a general scattering angle, 

m = / ( " /2 L S ; n \ - (3-33) 

' 1 + p(n/2) cos 2 6 ' 

just as for Thomson scattering. This formula is generally true for electric dipole transi- 
tions and so all we need compute is p(ir/2). 
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3.3.3. Doublets 

The next most complicated case is the Li-like doublet transition, eg CIVA1550. Here 
the ground level has J = 1/2 and there are two possible excited levels of which the lower 
energy level is J = 1/2. All transitions between sublevels are permitted and the net 
polarization is p(n/2) = 0. The higher energy excited level, associated with the shorter 
wavelength transition, has J = 3/2 and averaging over all of the permitted transitions 
gives p(w/2) = 3/7. If we average over both transitions according to their statistical 
weights, then we end up with p(n/2) = 3/11 Fig. (9). 

Lya, has a similar type of transition. Here the the wavelength separation of the two 
transitions is so small that the doublet will not be resolved and averaging over the two 
excited levels makes sense. However, in the case of CIV, the energy difference is larger 
than the likely thermal width, so that one could, for example, imagine continuum photons 
propagating out of an expanding flow, encountering the polarizing J = 1/2 — » 3/2 
transition first and then encountering the J = 1/2 — > 1/2 transition which erases all of 
this polarization. It is clear that the polarization is sensitive to the nature of BALQ 
outflows. 

3.3.4. Triplets 

The next simplest case is the B-like triplet transitions such as CII A1335. In this case, 
there are two choices for the ground state J = 1/2,3/2 and two for the excited state, 
J = 3/2, 5/2; the selection rules forbid direct transitions with J = 1/2 — ► 5/2. The 
energy difference between the two ground levels is small enough that they should be 
equally populated by collisions in the absence of radiative transitions. 

It is helpful to introduce the ionization parameter U which is the ratio of the num- 
ber density of hydrogen ionizing photons to the electron density, designated n. When 
the radiative excitation rate, ~ 10 4 £/nio s_1 ; exceeds the collisional excitation rate 
~ 600riio s , the ground state sublevels will be populated in an unequal fashion that 
must be computed by solving for all the transitions. (Note that both of these rates are 
likely to be much less than the spontaneous, de-excitation rate and so collisional de- 
excitation of the excited states is generally thought not to be an issue for permitted lines 
under AGN conditions.) Under these conditions, when the ground sublevels are radia- 
tively mixed, the resulting polarizations will differ from the values computed assuming 
as statistical population of the ground sub-levels due to collisions. In this case, the po- 
larization is increased from p(w/2) — 0.21 to p(ir/2) — 0.38. This increase in predicted 
polarization is typical. 

3.3.5. Supermultiplets 

The next level of complication arises when an ion in a single ground state can be excited 
into several different excited states under radiative mixing conditions. In order to solve 
for the population of the different sublevels and the polarization, we must consider these 
distinct multiplcts together. 

3.3.6. Magnetic Mixing 

A final variation, which is quite likely to be relevant in an AGN, arises when the 
magnetic field is strong enough that the cyclotron frequency, ujq = 1-8 x 10 7 (-B/1G) rad 
s _1 exceeds the radiative excitation rate. In this case, the relevant eigenstates are referred 
to the magnetic field direction rather than the normal to the scattering plane and density 
matrices have to be used to attack the problem in general. This is also known as the 
Hanle effect and is discussed at greater length here in the solar context by Dr. Stenflo. 
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An important consideration is the degree of Faraday polarization. The rotation of the 
plane of polarization is given by 

This is unlikely to be a factor at ultraviolet wavelengths, but could be significant in the 
optical. 

3.4. Interpretation 

3.4.1. Emission Line Clouds 

There is a standard model of the emission line gas based upon the notion of an emission 
line cloud, a stratified slab of gas of size ~ 10 13 cm, located at a radius R ~ 0.3 pc from 
the continuum source with a density ~ 10 10 cm -3 . The ionization state of the gas is 
determined by the relative importance of photoionization and recombination which is, 
in turn, controlled by the ionization parameter, U. Typically this is U ~ 0.1. These 
clouds have a photoionization temperature T ~ 10, 000 K and an equivalent sound speed 
~ 10 km s . However they are moving with speed ~ 10, 000 km s _1 and Mach number 
M = 1000 through a hotter and more tenuous confining medium. This is patently absurd! 
Nonetheless, this model does provide a good representation of the ratios of the observed 
line strengths. What is clearly required is a convincing dynamical model that retains 
the successful features of the atomic astrophysics. For the moment, we just consider the 
polarization in the context of the cloud model. 

If we consider a CIV A1550 photon propagating out of a "standard" cloud, the optical 
depth for solar abundance of carbon, mostly in a triply ionized state, is ~ 10 . Under 
these conditions, photons do not diffuse spatially out of the cloud as might, at first, be 
guessed. Instead, they undergo a random walk in frequency and escape when they migrate 
into the wings of the line where the cloud becomes transparent. The polarization really 
has to be computed using a Monte Carlo simulation and, under these conditions of high 
optical depth, it is not surprising that it is too small to be measured, even when there is 
velocity shear giving an anisotropic escape probability. The one conspicuous exception is 
the semi-forbidden (intercombination line violating the selection rule AS — 0) CIII] line 
A1909. Here the optical depth is closer to 10 and substantial linear polarization p ~ 0.01, 
was predicted and indeed has been reported in this line (Lee 1994, Cohen et al.1995, but 
see Ogle et aZ.1999). 

There is a rather different story when external photons are scattered by the emission 
line clouds. Here we expect a high albedo and roughly half the photons will undergo just 
one scattering (Korista & Ferland 1998). If the distribution of scatterers is anisotropic, 
then we might also expect to detect a linear polarization signal in permitted lines. This 
is not usually seen, which suggests that a particular line of sight contains at most one 
cloud at a given wavelength. This is a strong constraint upon models of the velocity 
distribution. 

3.4.2. Rayleigh Scattering by Lya 

Hydrogen is the most abundant element so scattering of Lya photons is likely to occur 
far into the wings of the line. The fine structure level splitting between 2P 3 / 2 and 1P\ii 
is quite small for this transition and off-resonance scattering is characterized by the 
(classical) Rayleigh scattering phase function (e.g. Stenflo 1980). This process is called 
Rayleigh scattering when the initial and final states of the atom or ion are identical. (This 
distinguishes it from Raman scattering which arises when they arc different.) Rayleigh 
scattering has been clearly observed in symbiotic stars (Nussbaumer, Schmid, & Vogel 
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1989). For large velocity shifts AT^, the scattering optical depth exceeds unity for 
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(3.35) 



(Lee & Blandford 2000). 

Of particular interest is the case of Lya, because in the damping wings the scattering 
phase function becomes that of the classical Raylcigh function, which enhances polariza- 
tion. The fundamental reason why this is the case is that near the line center, where the 
optical depth is large, most photons migrate in frequency space faster than they do in real 
space. They are therefore comparatively insensitive to the cloud shape and large scale 
velocity shear. However, in the wings of the line, the optical depth is much smaller and a 
large scale pattern in the cloud shapes, for example, translates into a measurable linear 
polarization. (Effects like this have been reported in high column density supershclls 
associated with starburst galaxies (Lee & Ahn 1998).) 

If the accretion disk of a quasar is warped so that some part of it is shadowed from the 
direct exposure to the central engine or has a thickness that increases slower than linearly 
with radius, a large column density Nhi may exist in the shaded region and Raylcigh 
reflection of Lya is expected. In this case we may expect up to 10 percent polarization 
in the Lya wings, which is consistent with ~ 7 percent polarization reported from the 
radio-quiet quasar PG 1630+377 (Koratkar et al.1995) 

3.4.3. Absorption Line Clouds 

A somewhat analogous situation is found for the absorbing clouds observed directly in 
the BALQs. These are believed to be located outside the emission line region where the 
ionization parameter U ~ 1 and the size is estimated to be even smaller than the size of 
the emission line clouds ~ (s/V) 2 R ~ 10 11 cm. (In a quite different type of model, it has 
been proposed that the emitting and absorbing gas originates from very much smaller 
radii and forms part of a space-filling flow, Murray et a/.1995. Many of the following 
considerations apply to this model as well.) 

The actual kinematics of line formation can be quite complicated. This is generally 
handled under the Sobolev approximation (Rybicki & Hummer 1978). An incident con- 
tinuum photon is scattered when it is resonant with a permitted transition taking place 
in the rest frame of the outflowing gas. The scattered line is redshifted in frequency by 
Vji/A. The surface occupied by gas resonant with a fixed observer frequency is called a 
Sobolev surface and can have a fairly convoluted shape. Some photons may encounter 
several Sobolev surfaces before escaping for good. The optical depth to absorption by 
ion X through a single Sobolev surface where the parallel velocity varies monotonically 
is 

dV\\ 
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t = I dsnxo = 
Substituting numerical values 



dV \\ 



ds 



(Is 



-l 
c 



v 



nx / dv aiv) . (3.36) 



T = 0.3X_ 4 f A a d y , (3.37) 
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where X = 10~ 4 A_4 is the abundance of the ion, / is the oscillator strength and 
dN2o/dVg is the hydrogen column density (in units of 10 20 cm~ 2 ) per unit velocity (in 
units of 10,000 km s _1 ). Resonance line scattering by the common ions occurs at col- 
umn densities about three orders of magnitude smaller than those required for electron 
scattering (Lee & Blandford 1997). 

We can compute the polarization expected from a particular model using Monte Carlo 
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simulations. It turns out to be possible to give factor 2 estimates for the polarization, 
by multiplying expressions which describe the most important factors in producing the 
integrated polarization, 

p ~ P (tt/2)D(t)AG , (3.38) 

In this equation, D(t) is a depolarizing factor that takes into account multiple scattering. 
It is typically roughly fit by an exponential D = exp(— r/b). A is a factor that takes into 
account the anisotropy in the escape probability. If the flow is uniformly expanding, it will 
be very difficult for the photons to escape in the radial direction and far easier for them to 
escape tangentially. This roughly doubles the polarization. Conversely, if the outflow is in 
the form of a jet, then small scattering angles will be favored with lower net polarization. 
Finally G is a geometrical factor which is supposed to account for the observed gas 
distribution. If the outflow is confined to an equatorial fan subtending a solid angle Af2, 
then G ~ Ail/3 and the electric vector will lie parallel to the projected symmetry axis. 
For a jet, the polarization will be perpendicular to the axis (which should be coincident 
with the radio axis). More extensive observations than have been possible so far will be 
needed to test the hypothesis that the lines are due to resonance scattering, through the 
dependence on atomic type (through p(7r/2)), and to decide upon the flow geometry. On 
this basis, the strongly polarized lines are expected to be HeIA584, OVA630, Ne IA736, 
NIVA765, CIIIA977, SiIIIA1206, A1IIIA1671 and MgIA2852. The modestly polarized lines 
include CIIA687, OIVA789 and NIIIA991 and weakly polarized lines include OVIA1034, 
HIA1216, NVA1240, SiIVA1396, CIVA1549 and MgIIA2798. CIIA858, NIIIA764, and 
OIVA609 are expected to be unpolarizcd and can be used to measure the amount of 
electron scattering. 

3.4.4. Polarization beyond the Lyman Edge 

Observations of the rest ultraviolet continuum from a few high redshift quasars have 
shown a strong polarization increasing irregularly to shorter wavelengths shortward of 
the Lyman continuum (Koratkar et aZ.1998). This may be as large as p ~ 0.2, although 
the observations were extremely difficult and are consequently a bit uncertain. One 
possible explanation is that there are several highly polarizing singlet lines in this region, 
like HeIA584, OVA630, NeIA736, NIVA765, or CIIIA977. If the outflow speeds associated 
with these relatively high ionization lines are large, V ~ 0.1c, then it is possible that 
the lines could overlap enough to give an apparent continuum polarization. It would be 
good to have the capability to repeat these observations. 

3.5. Summary 

• Spectropolarimctry provides a powerful diagnostic of the disposition of the broad 
emission and absorption line gas in quasars. 

• Resonance scattering should be variably polarized with the degree and direction 
dictated by fundamental considerations of atomic astrophysics and the flow geometry. 
The radiation observed in the troughs may be the scattered photons removed from other 
lines of sight. 

• By contrast, the optical continuum exhibits a fairly constant polarization, suggestive 
of electron scattering. 

• Emission lines are generally unpolarized, excepting the semi-forbidden line CIII] A1909. 
This is consistent with the standard cloud model of emission line formation. BAL troughs 
are variably polarized. 

• Large polarization rising with frequency has been reported to the blue of the Lyman 
continuum. This may be due to blends of strongly polarized, prominent singlet transitions 
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4. Neutron Stars 

4.1. Motivation 

Although physicists and astronomers (most famously Baade and Zwicky) were quick to 
appreciate the possibility that ~ 10 57 neutrons could assemble to form a self-gravitating 
neutron star, it was not until the discovery of radio pulsars in 1967, that there was 
compelling evidence that they really existed. To date we have cataloged over a thousand 
radio pulsars, know of hundreds of accretion-powered neutron stars in X-ray binaries, 
and are starting to find isolated neutron stars accreting from the interstellar medium. 
In addition, five radio, or rotation-powered, pulsars are observed to pulse at optical 
wavelengths (Chakrabarty & Kaspi 1998), at least seven as 7-ray pulsars (Thompson 
2000), and ~ 40 are detectable at X-ray energies (Becker, 2000). 

For the astronomer neutron stars are the most common result of evolution of a massive 
star. However, far from being an endpoint, they represent a rebirth often in a more 
luminous state than the progenitor star. For the physicist, neutron stars provide a 
magnificent cosmic laboratory, allowing us to witness the behavior of cold nuclear matter 
at supranuclear densities, the indirect effects of extremely high Tq conductivity and 
superfluidity, and, as we shall see, the consequences of magnetic field strengths perhaps 
nine orders of magnitude greater than we can sustain on earth. 

Neutron stars are also of special interest to the polarimctrist as they have already 
furnished the strongest and most rewardingly variable signals of any cosmic sources. It is 
possible to follow the change in the polarization (sometimes nearly completely polarized) 
through individual pulses from bright radio pulsars. It is also possible to study the 
average polarization properties of large samples of pulsars, viewed from a range of vantage 
points and, thereby build up a picture of the magnetic field geometry and try to determine 
the site of and the mechanism for their high brightness emission. Unfortunately, there 
is still no polarimetric capability at X-ray wavelengths where accreting neutron stars in 
binary systems emit most of their radiation. However, very strong linear polarization is 
anticipated and the details should be no less prescriptive of the emission. 

A third class of object, in addition to the accretion- and rotation-powered pulsars 
that is of particular interest at the moment is the magnetar. It appears that a minor- 
ity of neutron stars are formed with super-strong magnetic field ~ 10 14 — 10 15 G. As 
predicted by Thomson and Duncan (1995), these magnetars decelerate very quickly but 
still have a larger reservoir of magnetic energy that can be tapped to power 7-ray bursts. 
These field strengths are well in excess of the quantum electrodynamical critical field 
B c = m} e (? I eh = 4.4 x 10 13 G, where the cyclotron energy of an electron equals its rest 
mass. This, in principle, allows us to test the theory in a regime that is qualitatively 
quite different from that in which impressively high precision tests have already been 
made. (There is no real anxiety that the theory is suspect above the critical field but, 
as is the case with general relativity, there is a strong interest in performing the check.) 
Another class of X-ray source is the Anomalous X-ray Pulsars. These are possibly a late 
evolutionary phase of magnetars. 

In this section, we shall discuss the expected polarimetric properties of all three types 
of sources, although there are only observations of rotation-powered pulsars. 

4.2. Observation 

4.2.1. Rotation-Powered Pulsars 

Radio pulsars are spinning, magnetised neutron stars (e.g. Lyne & Smith 1998). The 
majority have spin period between 0.1s and 3s and surface magnetic field strengths 
~ 10 12 G, estimated from the rate at which they appear to slow down. The neutron stars 
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themselves appear to be mostly formed with masses quite close to the Chandrasckhar 
mass ~ 1.4 M©. Their poorly measured radii are ~ 10 km, consistent with there having 
central densities a few times nuclear as the best models of the nuclear equation of state 
imply. (We really do not know the interior composition at all well. It could be mostly 
neutrons or contain a large fraction of protons, hypcrons, pions or even free quarks. 
Accurate measurements of the radius along with the rate of cooling will provide important 
constraints on the equation of state of nuclear matter.) Other, impressive vital statistics 
of neutron stars include escape velocities ~ 0.3c, surface gravities ~ 10 14 cm s~ 2 and 
maximum spin frequencies (that are nearly attained in observed objects) ~ 1 kHz. The 
radio emission from pulsars has extremely high brightness temperatures which can, by 
some estimates, exceed ~ 10 30 K. 

The integrated pulse profiles of radio pulsars frequently show one, two or three pulses. 
This, and other observational evidence, has been interpreted in terms of an emission 
model where there is a strong "core" beam of emission close to the magnetic axis sur- 
rounded by a weaker "cone" beam. When the observer latitude is similar to that of the 
magnetic axis, a single, dominant core component is seen. Increasing (or decreasing) the 
observer latitude leads to a three peaked, cone-core-cone pattern. When the observer is 
more inclined to the magnetic axis, only the two cone components will be seen. These 
pulse profiles exhibit strong, broad band linear polarization varying through the pulse. 
Values p ~ 1 arc consistently measured at certain pulse phases in certain pulsars. The 
position angle swings regularly through the main pulse with a total swing that can be as 
high as ~ 180°. Strong circular polarization is also commonly measured near the center 
of the pulse, with the handedness often changing sign. 

Individual pulses, whose polarimetric properties can be measured in the strongest 
pulsars, are no less interesting. They show individual emission units, known as subpulscs, 
with durations typically a few degrees of pulsational phase. These can "drift" through 
the pulses appearing at progressively earlier or later phases in successive pulses. These 
subpulscs often exist in one of two orthogonal polarization states. Even shorter timescale 
features known as microstructure (or now even nanostructure) has been well documented 
and these too can exhibit high, though complex, polarization properties. 

The optical pulses from the Crab pulsar in the Crab Nebula have been particularly 
well studied. The pulse profile is cusp-like and the plane of polarization swings smoothly 
through ~ 70°, while the degree varies between p = 0.1 and 0.5. Another famous optical 
pulsar is associated with the Vela supernova remnant, shows similar strong, variable 
linear polarization. 

4.3. Physical Processes 
4.3.1. Curvature Radiation 

The large measured brightness temperatures imply that the radio emission is pro- 
duced by a coherent process. One widely discussed possibility is that the emission 
is some variant on coherent curvature emission whereby bunches of charged particles 
stream outward along the curving, roughly dipolar magnetic field lines from the star 
with ultrarelativistic speed (Lorentz factors 7 of several hundred) and radiate like giant 
electrons. The emission properties are like those already summarized for synchrotron 
radiation, with the important difference that the radius of curvature of the orbit, R, is 
energy-independent. The characteristic emission frequency is u> ~ 7 3 c/i?, lying in the 
radio band for R ~ 10 — 100 km. The radiation from an individual bunch is beamed 
within an angle ~ 7 -1 to the direction of motion. 

At a particular pulse phase, the observer sees emission from a curve though the mag- 
netosphere where the line of sight is tangent to the magnetic field. There is thought to be 
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a radius-to-frequency mapping so that the emission at a given frequency is concentrated 
over an interval of radius along this curve, and that this radius decreases with increasing 
frequency. The polarization from a tangent point will be quite strongly linearly polarized 
with electric vector parallel to the projected curvature vector on the sky. As the pulsar 
spins, this projected curvature vector will rotate on the sky and a characteristic swing of 
the plane of polarization will be produced. This is known as the rotating vector model. 
If we view an individual bunch from one side of its orbital plane then the other, we will 
see one sense of circular polarization followed by the opposite sense. This mechanism 
clearly has the ingredients to explain the radio polarization observations. However, it is 
a bit puzzling how a totally polarized pulse can be formed in this manner. One possi- 
ble explanation is that the bunches are quite strongly flattened and they radiate most 
strongly perpendicular to their flattening plane, where the polarization will be most 
strongly linear. 

4.3.2. Maser Processes 

There has recently been a resurgence of interest in maser emission models. It is rela- 
tively easy to imagine that the necessary population inversion will develop in the outflow- 
ing plasma. A typical pulsar can develop an EMF of ~ 10 13 — 10 16 V and a small fraction 
of this potential difference developing in a transient "gap" will create a fast stream of 
electrons and/or positrons that can stream through more slowly moving (though still 
ultrarelativistic) plasma. What is a bit harder is to find suitable wave-particle interac- 
tions that can lead to an overall negative absorption coefficient. As an illustration of the 
difficulty, note that synchrotron radiation in vacuo has an absorption coefficient 

' / "dA^^- 7 V(^7)- (4-39) 
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The single particle emissivity, p v (v,"/) is proportional to ^ 1 / 3 -y -2 / 3 at low frequency and 
increases with energy at high frequency. Hence, synchrotron absorption is necessarily 
positive and maser action is precluded, independent of the particle distribution function. 
Similar conclusions have been drawn for other emission mechanisms in vacuum. 

However this conclusion does not necessarily follow if there is a plasma present. One 
particularly interesting case is the so-called anomalous cyclotron (otherwise known as 
cyclotron-Cerenkov) resonance between a wave with angular frequency 10 and wave vector 
k interacting with an electron (or positron) moving with velocity v, in a field where the 
non-relativistic gyro frequency is wg. This occurs if 

uo - k\\V\\ = -ujgM ■ (4.40) 

where || refers to the component along the magnetostatic field. 

This equation needs some interpretation (Lyutikov, Blandford & Machabeli 1999). 
Consider, for simplicity, a circular polarized wave propagating along the magnetic field. 



A resonance satisfying the condition Eq. (4.40) clearly requires that the phase velocity 



of the wave be less than c. If we transform into the guiding-center frame of the electron, 
then the wave angular frequency changes sign indicating that it is propagating in the 
opposite direction along the magnetic field, with the same sense of circular polarization. 
This means that it resonates with particles gyrating in the opposite sense around the 
field (i.e. with opposite charge) than is the case in a regular cyclotron resonance. Put 
another way the particle outruns the wave so that, in the frame where the wave is at 
rest, the particle follows the electric vector as it spirals around the magnetic field. A 
consequence is that, in the electron guiding-center frame, a quantum of wave energy has 
negative energy. Therefore in exciting the electron to a higher state of gyration it emits 
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Figure 10. Location of core and conal emission regions in a maser model of pulsar emission 

(after Lyutikov et al. 1999). 



a quantum and vice versa. (In the rest frame, population inversion now requires having 
more particles in the lower gyrational state!) Not surprisingly, this arrangement can lead 
to maser action in the outer magnetosphere and this has been proposed as a mechanism 
to produce the core emission. The modes arc naturally circular polarized, although the 
handedness depends upon the details of the electron and positron distribution functions. 
This is not the only possible way to have a maser process. There is a second resonance 
associated with the curvature drift motion ^drift of the gyrating electron as it moves along 
the curving magnetic field 

LO - k\\V\\ = fcj_Vdrift. (4-41) 

This will be perpendicular to the curvature plane and will consequently produce emis- 
sion with polarization orthogonal to that predicted by the rotating vector model. This 
mechanism has been invoked to account for the cone emission (Fig. (10)). 

4.3.3. Propagation Effects 

There is unfortunately a complication (Arons & Barnard 1986, Lyutikov et al.1999, 
Hirano & Gwinn 2001). The emitted radiation must propagate through the outer mag- 
netosphere. This can lead to genuine absorption at the normal cyclotron resonance by 
more slowly moving electrons. Landau damping is also a possibility. Furthermore, there 
are refractive effects that may imprint additional polarization on the emergent radiation 
in much the same way as occurs in the ionosphere. There can be mode conversion, for 
example from a subluminal ordinary mode to a propagating electromagnetic wave. Fi- 
nally, and perhaps most interestingly from a physics perspective, there can be non-linear 
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scattering effects (Lyutikov, 1998). At the high brightness temperature (or, equivalently, 
large occupancies of individual quantum mechanical states) found in pulsar radiation 
there will be a large amount of scattering between different radiation beams. These can 
be mediated by individual electrons, in which case the interaction is known as induced 
Compton scattering, or by collective wave modes of the plasma (known as Raman scat- 
tering when the scatterer is an electrostatic wave). There is not space to discuss these 
rather complex processes further, save to remark that, the associated matrix elements 
have a strong sensitivity to polarization and frequency and that this can be used to 
identify them. The theory is starting to match the observations in its richness! 
Some of these propagation effects can also be relevant in the interstellar medium. 

4.3.4. Thomson Scattering in a Strong Magnetic Field 

The cross section for Thomson scattering must be changed if there is a strong magnetic 
field present (e.g. Meszaros 1992). In the limit when the wave angular frequency u> << 
u>g, or E « 44(i?/10 12 G) keV, the electrons are constrained to move along the magnetic 
field like beads on a wire. The dominant cross section is between polarization states with 
k,E,B coplanar. It is clearly given by 

^- = r 2 sin 2 <9sin 2 6>', (4.42) 

ail 

where 0, 8' are the angles made by the incident and scattered wave vectors with the 
magnetostatic field. 

4.3.5. Inverse Compton Scattering in a Strong Magnetic Field 

This anisotropy in the cross section introduces an additional complication to inverse 
Compton radiation. This is because the incident photon propagates in a direction mak- 
ing an angle ~ 7 -1 to the magnetic field and so the total scattering cross section as 



given by Eq. 1.42 is reduced by a factor ~ 7 2 from the Thomson value assuming that 
the frequency in the electron rest frame (a;') is less than loq. However, under these 
circumstances, we must also consider the effect of the "E x B" drift of the electron per- 
pendicular to the magnetic field. This will produce an oscillatory electron motion that 
is larger than the field parallel motion by a factor ju>'/u>g- The azimuth and incident 
polarization-averaged cross section will be given by 

^ = !(^) 2 (l + cos 2 #)(l + cos 2 tf'), (4.43) 

for incident frequency oj satisfying ioqJ^uj < loq. In practice, this leads to a rather 
complex polarization pattern. 

4.3.6. Quantum Electrodynamical Effects 

The virtual electron-positron plasma that comprises the QED vacuum affects the prop- 
agation of radiation through it. These are clearly likely to be of importance when the 
magnetic field strength approaches the critical field strength, B c . 

An external magnetic field causes photons of different polarizations to travel at slightly 
different speeds. The CP-invariance of electrodynamics tells us that the two modes must 
be linearly polarized. Specifically, the indices of refraction of both modes differ from unity 
and are given by (Heyl & Hernquist 1997a): 

a 8 /5±\ 2 

n - = 1+ i^45 U: + -« (444) 



n|| = l + 73^sin 2 6»^ + --- . (4.47) 
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7t.ii = 1 H — H , (4.45) 
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for hco <C m e c 2 and B <C B c . B± is the component of the magnetic field perpendicular 
to the propagation direction of the photon. A photon in the perpendicular polarization 
has its electric field vector perpendicular to the projection of the magnetic field into the 
transverse plane, and similarly for the parallel polarization. 

For fields stronger than B c , the index of refraction for the mode with the electric field 
perpendicular to the external magnetic field saturates at 

n± = 1 + — -sin 2 0-1 , (4.46) 
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while the index for the other mode increases without limit 
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Therefore, there are two natural limits to the behavior. In the weak field regime the 
vacuum polarization may be sufficiently strong to decouple the polarization states as 
they propagate through the magnetosphcre (Hcyl & Shaviv 2000, 2001). Normal pulsars 
as well as strongly magnetized white dwarfs fall in this regime. In the strong-field regime 
magnification and distortion of the image of the neutron star surface may become impor- 
tant in addition to the decoupling of the propagating modes (Shaviv, Heyl & Lithwick 
1999). 

The processes of photon splitting and one-photon pair production are forbidden in field- 
free regions, the first by Furry's theorem and the second by four-momentum conservation. 
However, in the strong magnetic fields surrounding a neutron star, both processes may 
be important. Photon splitting most strongly affects photons in the _L-mode which may 
split into two photons in the ||-mode (Adler 1971, the mode-naming convention used here 
is opposite to that used by Adler); this both distorts and polarizes the photon spectrum. 

One-photon pair production has an energy threshold of Huj sin > 2m e c 2 for photons 
in the ||-mode. The threshold for photons in the _L-mode is slightly larger TiuismO > 
m e c 2 (l + y/l + 2B/B C ). Near the thresholds, in strong magnetic field especially, the 
cross-section for this process is complicated by the formation of the pair in discrete 
Landau levels or a positronium bound state (Daugherty & Harding 1983; Usov & Melrose 
1996). 

4.4. Interpretation 
4.4.1. Rotation-powered Pulsars 

Partly because it provides such a straightforward interpretation of the polarization 
data, the curvature radiation model is probably the favorite explanation for pulsar radio 
emission. However, it does have some drawbacks. One of these is that it is difficult to 
maintain a compact bunch for very long as the electrons travel along trajectories with 
different radii of curvature and are subject to radiation reaction. Another problem is that 
it has proven hard to find a suitable plasma instability which will allow charge particle 
bunches to grow. A third and currently controversial observational claim is that the 
emitting area is much larger than expected if the bunches form in the inner magnetosphcre 
as the rotating vector model requires. Specifically, Gwinn, et aZ.(2000) find that the size 
of the Vela pulsar emission region is roughly ten per cent of its light cylinder radius, 
~ c/Cl, suggesting that the emission originates in the outer magnetosphere. Conversely, 
Cordes (2001) finds that the source is unresolved. 

The maser explanation, by contrast, only works in the outer magnetosphere and the 
average direction of polarization should be orthogonal to the projection of the pulsar 
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spin axis on the sky. This can be tested using the Crab and Vela pulsars where X-ray 
jets are observed which, although their formation is not understood, are presumed to be 
along the projected spin axis. In the case of the Crab pulsar, the situation is ambiguous 
because it is unclear if the two pulses come from one or two magnetic poles. However, 
in the case of the Vela pulsar, the electric vector is unambiguously perpendicular to 
the projected spin axis, consistent with the maser model. It could also be consistent 
with curvature radiation if propagation effects are important. It should be possible to 
discriminate between these two models observationally. 

The pulsed optical radiation and X-ray radiation seen from several radio pulsars is 
generally thought to be incoherent synchrotron and inverse Compton radiation, respec- 
tively. However, the location of the emission region, and, in the case of the X-rays, the 
source of the incident photons (the surface of the star or coherent radio emission from 
the magnetosphere) is currently undecided. Suffice it to remark it here that polarization 
arguments figure prominently in these debates. 

4.4.2. Accretion-powered Pulsars 



In addition to being strongly polarizing, the opacity, Eq. |4.42| , is highly anisotropic. 
So, even though we cannot measure the X-ray polarization directly, at present, it does 
have a strong effect on what we observe. In particular, the mass that accretes onto a 
spinning neutron star with a surface field ~ 10 12 G is likely to be channeled toward 
the magnetic poles and, if the accretion rate is large enough, there will be a significant 
Thomson opacity at the poles. However, it will be much easier for the radiation to escape 
along the direction of the magnetic field than in the transverse direction. For this reason, 
X-rays are thought to emerge in two broad pencil beams about the magnetic axis as is 
observed. 

4.4.3. Surface Emission from Isolated Neutron Stars and Magnetars 

The atmospheres of neutron stars are thought to emit strongly polarized radiation 
(Pavlov & Shibanov 1978). The opacities in the two polarization modes of the atmo- 
spheric plasma may differ by several orders of magnitude (Lodenqual et al.1974). The 
opacity in the extraordinary mode (i.e. E _L B) is generally a factor of (w/wg) 2 smaller 
than in the ordinary mode. Since the atmospheres are typically at a temperature of 
several million degrees, the natural place to study this emission is in the X-rays. Fur- 
thermore, as we shall see, the vacuum significantly affects the propagation of radiation 
passing through it at X-ray and higher energies. 

Although the emission at the surface may nearly be fully polarized, one observes radi- 
ation from regions with various magnetic field directions. In this vein, Pavlov & Zavlin 
(2000) argue that the net polarization in the X-rays is on the order of ten percent and 
decreases for more compact stars. However, this treatment ignores the fact that QED 
renders the vacuum birefringent. The field strength varies sufficiently gradually, that is 

W <<|at| ' ( " 8) 

where 

mss ±L?L(*£\\ (4.49) 
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in the weak field regime, that the two polarization modes are decoupled. Radiation 
produced at the surface with its polarization direction perpendicular to the local magnetic 
field direction will keep its polarization perpendicular to the field even as it passed 
through regions where the field direction changes. The observed polarization reflects 
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the direction of the field at a distance 

/ \ 2/5 -i /c- 

r«1.2xl0 7 ^ A (— £ ) * (sin/?) 2/5 cm, (4.50) 

Vl0 30 GcmV Vl0 17 HJ V l! ' V ' 

from the center of the star (Heyl & Shaviv 2000, 2001). Here /i is the magnetic dipolc 
moment of the neutron star, and (3 is the angle between the dipole axis and the line of 
sight. QED ensures that the strongly polarized radiation at the surface of the neutron 
star remains strongly polarized until it is detected; therefore, the simple detection of 
strongly polarized X-rays from the atmosphere of a neutron star will verify a thus far 
untested prediction of QED. Cheng & Ruderman (1979) used a similar argument to 
account for the strong polarization of radio emission from pulsars. 

Although, the QED process of one-photon pair production plays a crucial role in radio 
pulsars by fueling the plasma that produces the emission (e.g. Daughcrty & Harding, 
1982), the threshold for the reaction is much higher than the typical energies from the 
surface emission. The cross-section for photon splitting increases dramatically with in- 
creasing photon energy, oc E 6 (Heyl & Hernquist 1997b) and is only important above 
10 keV even in the strongest magnetized sources. 

4.5. Soft-Gamma Repeaters 

Thompson and Duncan first argued that the soft-gamma repeaters are neutron stars 
fueled by a dynamic magnetic field whose strength greatly exceeds B c . In their quiescent 
state, these objects emit thermal radiation from their surfaces and the discussion of the 
previous subsection applies. If their surface fields are sufficiently strong (the surface 
field is expected to exceed the value estimated by spin down of ~ 10 15 G), magnetic 
lensing may be important for photons whose polarization is parallel to the magnetic field 
(i.e. the ordinary mode); however, thermal emission in this mode appears to be strongly 
suppressed. 

However, what makes the soft-gamma repeaters unique is that they burst. In fact, the 
soft-gamma repeater, SGR 1900+14, is the only object beyond our solar system to have 
had contemporary geophysical consequences (it ionized the nightside upper atmosphere 
nearly to daytime levels). This soft gamma-ray emission is generally well below the 
threshold for one-photon pair production, but photon splitting should degrade the ener- 
gies of the photons by at least a factor of two and polarize them by converting photons 
in the extraordinary mode to the ordinary mode (Baring & Harding 1997). Observing 
this tracer of photon splitting would require gamma-ray polarimetry. 

4.6. Summary 

• Magnetized neutron stars provide cosmic laboratories where we can observe unique 
polarization effects in action and use them to identify the emission mechanism. 

• Radio pulsars offer the richest polarization data set outside the solar system. They 
are strongly diagnostic of the emission mechanism and the effects of propagation. 

• Accretion-powered pulsars introduce new effects associated with strong field anisotropic 
emission and scattering in the ~ 10 12 G surface fields. Even though the X-ray polariza- 
tion is not yet measured, it is important in determining the total spectrum and pulse 
profile. 

• X-ray pulsars allow us to address important physics questions, like the composition 
and compressibility of cold matter at supra-nuclear density. 

• The simple detection of strongly polarized X-rays from the atmosphere of a neutron 
star (or optical radiation from neutron stars with B > 10 13 G) will verify the prediction 
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that QED renders the vacuum bircfringent and provide an estimate of the radius of the 
star itself. 

• The convincing case that magnetars exist with surface fields well in excess of the 
critical field, (4.4 x 10 13 G), offers the equally exciting, (though observationally very 
challenging) prospect of testing quantum electrodynamics in a regime far removed from 
terrestrial investigation. 



5. Black Holes 

5.1. Motivation 

There is now very good evidence for the existence of black holes in the universe. They 
appear to be a common endpoint of the evolution of massive stars in our Galaxy and 
nearby galaxies and we know of roughly ten good cases where the dynamically determined 
mass significantly exceeds the Oppenheimer-Volkoff limit (for neutron stars) ~ 2.5 M© 
and the Chandrasekhar limit (for white dwarfs) ~ 1.4 M©. In addition, there are many 
more cases of transient X-ray sources where the circumstantial evidence, in the absence 
of dynamics, is pretty convincing. A significant fraction of massive stars must end their 
life this way. 

Similarly, dynamical studies of the nuclei of nearby galaxies reveal the presence of 
"massive dark objects" which, if they were, for example, clusters of compact objects, 
would be very short-lived. Identifying them with massive black holes is by far the most 
conservative conclusion to draw. It appears that the nuclei of most normal galaxies, 
including our own, contain black holes with masses in the range ~3x 10 6 — 3 x 10 9 M©. 
(There are speculative suggestions that there may be a large population of intermediate 
mass black holes, perhaps relics of the first generation of stars.) 

The existence of black holes is, arguably, the most far-reaching implication of the gen- 
eral theory of relativity. The theory has been probed in the weak field regime and passed 
all quantitative tests with an accuracy that can be as small as ~ 3 x 10 -4 . It is in 
the nature of the theory that, if we understand the laws of physics under these circum- 
stances, it is simply a question of geometry to describe strong field environments, when 
the equivalent Newtonian potential approaches c 2 . If the theory in its essential simplicity 
is correct, then the metric of an asymptotically flat black hole spacetime (excluding some 
mathematical niceties) is essentially given. Indeed, in one of the greatest successes of 
mathematical physics, we have a closed form version of the metric of a spinning black 
hole, known as the Kerr metric, and essentially all classical physics that can be discussed 
in a flat spacetime can also be discussed around a black hole; there are no difficulties of 
principle. There are, however, considerable difficulties in execution (which have mostly 
been overcome using numerical calculations.) Although we know of more general space- 
times that include a gravitationally significant charge or orbiting mass, we believe that 
these are irrelevant to observed black holes and that astronomers need only be concerned 
with the Kerr metric. 

However, it is logically possible that the theory of general relativity could be wrong 
or incomplete and that, as a consequence, black holes are fundamentally and observably 
different from their general rclativistic description. For this reason, it is vitally important 
that we try to find ways to probe the spacetime around black holes, now that we know 
where to find them. In this regard, observing black holes provides a far more telling 
test of relativity theory than cosmological observations. This is because cosmological 
observations are seriously compromised by our deep ignorance of the nature of dark 
matter and energy as well as the effects of evolution. 
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We already know that gravitational waves exist. Binary pulsars are observed to lose 
orbital energy at rates that agree with theory to a fraction of a percent. However, this 
mostly tests linearized theory even in the sources. The direct detection of gravitational 
waves, which we hope will happen one day, is also is only a linear test. The ultimate 
test of general relativity is to make detailed observations of gravitational waves from 
coalescing black holes; an observation that I suspect will be not be technologically fea- 
sible for some time. In addition to testing strong field relativity, this can also provide 
much useful astrophysical information on galaxy merger rates, AGN evolution and so on. 
Computing the wave forms in necessary generality is a major challenge to computational 
science. What is relevant in the present context is that gravitational waves have natural 
polarization states, just like electromagnetic waves, and much of the information from 
these coalescences will be encoded in the polarization details. (There are other strong 
field sources of gravitational radiation that have been considered, notably topological 
defects like cosmic strings. Unlike the case with black holes, there is no observational 
evidence yet for their existence. However, if they are ever discovered, then it may well be 
their gravitational radiation polarization that is their distinctive signature. Computing 
this polarization is a good project which appears to have been mostly ignored.) 

There is a second and quite different reason for being interested in black holes. This 
is that we do not understand properly how they work. We have already introduced jets, 
disks, and outflows; and discussed how polarization observations can teach us about their 
properties. We already know, from direct observation, that all three of these continue 
down to relatively close (in logarithmic terms) to the black hole. This is where most of 
the energy is released. However, we do not understand how this all happens and how 
the flow around the black hole depends upon the mass and the spin of the hole as well as 
the accretion rate and the immediate environment. It appears that the answers to these 
questions will only be found by exploring the black hole itself. 

A third contextual aspect of this study is that we are beginning to suspect that black 
holes have a much larger and more active role in galactic and extragalactic astronomy 
than used to be the case. It is increasingly likely that gamma ray bursts are associated 
with the formation or augmentation of black holes and that these have major environ- 
mental impacts on their surroundings and could soon become useful cosmological probes. 
Black hole transients provide the dominant hard X-ray emission of galaxies like our own 
and create powerful outflows. The discovery of dormant, or near-dormant black holes 
in the nuclei of normal galaxies has affirmed the long-standing black hole model of ac- 
tive galactic nuclei, including quasars and giant, double radio sources. However, the 
implications have much broader implications than the properties of AGN per se. The 
quasars themselves provide the best cosmologically distant beacons that we have and 
they allow us to study the intergalactic medium, matter and cosmography. Furthermore, 
it is becoming increasingly apparent that they have an much more active role in the very 
formation of galaxies, both in the initiation and perhaps in the cessation of the process. 

The long term observational goal, then, is to verify that black holes are described by 
the Kerr metric and to measure their masses and spins in such a way as to elucidate their 
role in stellar and galactic evolution. In this section, we will try to show how polarization 
observations can contribute to meeting this objective. 

5.2. Observation 

We have already described most of the relevant observations of black holes including the 
fairly strong dynamical measurements of their masses and the Fe K line emission which 
provides the strongest evidence to date that black holes spin relatively rapidly. The 
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most relevant, existing observation for the purpose of this lecture are the measurements 
of linear and circular polarization. 

5.2.1. Sgr A* 

The center of our Galaxy appears to be identified with the radio source Sgr A* . There 
is now excellent dynamical evidence that is a "dark, compact object" with a mass 2.6 x 
10 6 Mq and a black hole is by far the most conservative interpretation. The source 
is nearly at rest and stars can be tracked moving (and accelerating) around it (Ghez 
et aL2000). The source has a spectrum which peaks at ~ 300 GHz and the image 
is broadened at radio wavelengths by interstellar scattering. At 43 GHz, the average 
scatter-broadened size is reported to be ~ 3 x 10 13 cm, (Lo et aL1999) and to scale 
roughly oc A 2 . Sgr A* is a weak and soft X-ray source. Interestingly, a 106 d periodicity 
in the radio emission has also been reported (Zhao, Bower & Goss 2001). 

Sgr A* has long been known to have negligible linear polarization at radio and mm 
wavelengths. This is not a surprise because the Faraday rotation is expected to be quite 
high so that the differential (in both angle and frequency) rotation is also large enough 
to depolarize all measurements. What is a surprise (as discussed here by Hildebrand) is 
that Aitken et aZ.(2000) measure ~ 10 percent linear polarization at 150 GHz, using the 
SCUBA instrument on the James Clerk Maxwell Telescope, although there is an upper 
limit of ~ 1 percent at 86 GHz (Bower et aZ.2000). Clearly there is a need to confirm 
the SCUBA measurement. 

More recently, it has been discovered that Sgr A* exhibits quite strong circular po- 
larization, (Bower, Falcke & Backer 2000). The 5 GHz degree of circular polarization 
appears to have been stable at a value of ~ —0.003 for nearly twenty years. At higher 
frequencies, up to ~ 43 GHz, the degree of circular polarization appears to increase up 
to a few percent and become increasingly variable, doubling in a few days. 

5.3. Physical Processes 

5.3.1. Spinning Black Holes 

The spacetime around a spinning black hole is described by the Kerr metric expressed 
in Boyer-Lindquist coordinates with G = c = 1 

ds 2 = - (1 - 2mr/p 2 )dt 2 - (4amr sin 2 9/p 2 )dtd<j) + (p 2 /A)dr 2 + p 2 d9 2 

+ (r 2 + a 2 + 2mra 2 sin 2 9 / p 2 ) sin 2 9d<j) 2 , (5.51) 

where 

p 2 = r 2 + a 2 cos 2 9 , (5.52) 

A = r 2 - 2mr + a 2 , (5.53) 

and m is the mass such as would be measured by the orbit of a distant satellite, and 
a < m is the specific angular momentum of the hole, as could be measured operationally 
by the precession rate of a gyroscope. 

There is an event horizon, 7Y, which is located where A = i.e. where the radial 
coordinate r = r + = m + (in 2 — a 2 )z. Particles on timelike or null geodesies must be 
inwardly moving within r + which leads to the interpretation that H. represents a surface 
of no return. The four velocity, u — {dt/dr, dr/dr, d9 / dr, dfi/dr} of a material particle 
satisfies 

g a0 u a u fi = -l. (5.54) 

The equation of a photon, following a null geodesic is given by g a .j3dx a dx 13 = 0, supple- 
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merited with equations representing the conservation of energy and angular momentum 
as well as an additional integral of the motion. 

The angular velocity SI — d<j)/dt of a particle, orbiting with hxed r, 8, therefore satisfies 



u 



.goo + 2S1.9O0 + ft 94,4,} = -1 • (5.55) 

This implies that fi m [ n < ft < ft m ax where f2 m ; n > when r + < r < r e = m + (to 2 — 
a 2 cos 2 6) 1 ' 2 . The radius r e is known as the static limit and the region between it and 
the horizon, where inertial frames are dragged by the spin of the hole, is known as the 
ergosphere. A particular significance of the ergosphere is that orbits of negative energy 
(including rest mass) exist within it. As r — ► r + at the event horizon, 

ftmi„ -> ftmax -► Off = a/(r 2 + + a 2 ) 1 ' 2 , (5.56) 

the angular velocity of the hole. 

A remarkable theorem due to Hawking states that the area of the horizon, which 
can be computed from the metric to be A — J-yAgeeg^) 1 ^ 2 'd9dcp = 47r(r"j_ + a 2 ) cannot 
decrease. We can use this to define a so-called irreducible radius r^ and irreducible mass 
mo through 

r = 2to = (A/4tt) 1/2 . (5.57) 

This immediately implies that a = r 2 iln- It turns out that the area is proportional to the 
thermodynamic entropy. Now imagine that we exchange some mass and some angular 
momentum with the hole, reversibly (and therefore at constant area) from just outside 
the horizon. These must be added according to 

dm = Cld(am) . (5.58) 

This equation can be integrated to give 

m =[l-(ft,o) 2 P/ 2 ' (5 ' 59) 

Imposing the condition a < to, we find that there is a mass to — too < 0.29m, that 
can, in principle, be extracted from the hole. The main way that this is thought to 
occur naturally is through the agency of large scale magnetic field that threads the event 
horizon of the black hole. This magnetic field can exert a torque on the hole, similar to 
the magnetic torques acting upon the sun and neutron stars, for example. Black hole spin 
provides a plausible power source for high energy phenomena like ultrarclativistic jets 
and gamma ray bursts and this is one reason why black holes are commonly thought to 
be spinning rapidly. (Even if the spin is not a significant power source, then the specific 
angular momentum of the gas that accretes onto a black hole is generally so large that 
it is very hard to imagine slowly spinning holes ever being formed.) 

For present purposes, though, what is most important is that, in a rapidly spinning 
hole, the accreting matter can form a disk extending quite close to the horizon. Specifi- 
cally, if we consider circular Keplerian orbits around a hole then these are stable down to 
a radius of marginal stability which is located at 6m. for a non-rotating (Schwarzschild) 
hole and approaches the horizon as a — > to. In addition, it is possible for strong pressure 
gradients within the disk to support matter in non-Keplerian orbits inside 6to. This 
means that gas may survive quite a long while in and around the ergosphere before 
crossing the horizon or being ejected. The dominant emission may come from this re- 
gion, and as the radiation escapes, its trajectory and the propagation of its polarization 
can be significantly influenced by the curvature of the spacctimc. This provides us with 
a potential probe of the Kerr metric. 
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5.3.2. Geometrical Optics of Plasma Waves in Flat Space 

We are interested in the propagation of plasma waves in the curved spacetime around 
a black hole. However, for the moment, let us consider the propagation of waves in 
flat space under geometrical optics. This is appropriate because the wavelengths that 
we are considering, at least for electromagnetic radiation, are always much smaller than 
the horizon radius. It is convenient to exploit the analogy with Hamiltonian particle 
dynamics. Under the eikonal approximation, we can define a phase <f> such that V(/> = k 
and d(f>/dt = — uj. The phase velocity is, as usual, defined by Vj, = to/k. We assume the 
existence of a dispersion relation 

u) = tt(k,x,t). (5.60) 

Equivalently, there is a Hamilton-Jacobi equation 

^+ft(V0,z,i)=O, (5.61) 

which must be satisfied. The three Hamilton equations are 

dk 
lit 
dx duj 

dt ~ dk 
dui d£l 
~dt ~ ~dt 
where 

d d 



Vil, (5.62) 

= V g , (5.63) 

(5.64) 



* 7 flf^- V ' (5 ' 65) 

and V g is recognized as the group velocity. These three equations govern the propagation 
of plasma modes in a spatially inhomogencous and temporally varying medium, under the 
short wavelength approximation. In a cold, unmagnetized plasma the dispersion relation 
is fl = (bjp + ^k 2 ) 1 / 2 and V g — c 2 /V<f,. We can think of wave quanta - plasmons - and the 
energy they carry, as moving along a path x(t) with the group velocity. These plasmons 
are conserved; i.e. the wave energy density U can be shown to obey a conservation 
equation of the form 

Consider, for example, shear Alfven waves. The dispersion relation is u> = k-B /(Airp) 1 ' 2 
and the group velocity is V g = B/^np) 1 / 2 . The wave packets propagate along the mag- 
netic field along with the energy although k can be directed at a large angle to B. 

The propagation of the polarization can be most simply approached by decomposing 
the given wave into its normal modes, propagating each along the direction of the group 
velocity, and compute the relative change in phase (to lowest order in the eikonal approx- 
imation this is simply J dx ■ k along the path.) Of course the character (i.e. polarization, 
local phase velocity etc. ) of these modes will change, but provided we are in the WKB 
limit, the modes are distinguished, non-degenerate (see below), and there is no mode 
crossing (which can occur and has to be handled more carefully), these changes will 
change adiabatically. As a result, they can be tracked and the total phase difference 
along a path can be computed. 

The total flux can be computed by using the conservation of intensity along the path. 
Equivalently, we say that the phase space density of individual quanta of wave excitation, 
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in individual modes is conserved along paths. If there is emission or absorption along 
the path then it is straightforward to write down the equation of radiative transfer and 
use the local emission and absorption coefficients to evolve the intensity (e.g. Rybicki & 
Lightman 1979; Bekcfi 1966). 

In practice, of course, all of this can easily become quite involved. However it is 
important to understand the principles because these alert us to the sort of effects we 
might expect to observe. 

5.3.3. Magnetized Accretion Disk 

As a more pertinent illustration of some of these ideas, let us consider electromag- 
netic wave modes propagating through an accretion disk containing a strong, though 
disordered, magnetic field. 

Consider a magnetoactive plasma with X = ujp/LU 2 ,Y — ujgJlo < 1 where u>p is 
the plasma frequency and log is the electron gyro frequency. Under so-called "quasi- 
longitudinal" conditions - essentially when cos6> < Y, where 9 is the angle between k 
and B, the electromagnetic eigenmodes are elliptically polarized with axis ratio r = 
1 ± Y sin 9 tan 9 and phase velocity difference AV = cXY cos 9. 

Now suppose that synchrotron (or cyclotron) radiation is emitted within an accretion 
disk of thickness H. The major axis of the polarization ellipse will be Faraday rotated at a 
rate A$/ds = AVlu/2c 2 . Now if, as we expect, the magnetic field direction reverses often 
along a ray and if, as also anticipated, X\Y\ioH/c >> 1, then we expect that the emergent 
linear polarization will be vanishingly small (cf §1.3.4). (The limiting polarization along 
an individual ray is likely to be determined by the decrease in the density rather than 
the magnetic field strength.) 

The circular polarization is a bit more problematic. If we suppose that there is a net 
magnetic field normal to the disk, as is true of some models, then there should be a 
preferred sense of the circular polarization that is emitted in cyclotron or low energy 
synchrotron radiation. This will be largely preserved in propagating out of the disk. 

If the emitted radiation is effectively unpolarized, we can analyze the production of 
circular polarization due to Faraday conversion by decomposing each wave into the eigen- 
modes. Thus consider a single eigenmode propagating out of the disk through a spatially 
varying magnetic field. We suppose that the variation happens relatively slowly on the 
scale of the wavelength so that the polarization ellipse adjusts adiabatically (with no 
mode crossings). Next, suppose that the two eigenmodes are launched with equal ampli- 
tude and that the field is uniform. The beating between the two eigenmodes will result in 
a circular polarization of amplitude Y sin 9 tan 9 that changes sign as the plane of linear 
polarization rotates. If cither the Faraday depth is large, or the sign of the magnetic 
field is as likely to be negative as positive, then the limiting circular polarization that 
emerges from the disk is equally likely to have either sign and so there will be no net 
circular polarization. 

Now let the magnetic field direction vary along a ray. If the angle 9 varies, then there 
will be a corresponding change in the axis ratio of the polarization ellipse, but still no 
preference for one sign over the other. However, if the azimuthal angle <p relating B 
to k changes in a systematic fashion along all rays then a net phase difference between 
the two modes will develop. In the same way that Faraday conversion creates circular 
polarization, such a phase difference will also create circular polarization. The main 
difference is that while Faraday conversion depends upon the direction of the magnetic 
field, and hence will not lead to a net polarization for randomized fields, the new mech- 
anism depends upon the rate of shearing and only the strength of the field. Thus, it is 
possible to conceive of situations in which the rate of shearing and the typical length 
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scales over which the magnetic field reverses are related in such a manner that a net 
circular polarization is produced without a commensurate linear polarization. 

This situation is precisely what might be anticipated in a magnetized accretion disk. In 
the disk interior, the typical field direction will trail to reflect the differential rotation in 
the disk. However, the field will be swept back by progressively smaller angles as the ray 
approaches the disk surface, corresponding to a net rotation of the average azimuthal 
angle <f>. There will only be a preferred sense of limiting circular polarization, if the 
magnetostatic field is still changing in this systematic manner over the last radian of 
Faraday rotation. The net circular polarization will be ~ clnY/ouHX. 

These are some of the subtle effects that could be present in an accretion disk and which 
could, under some circumstances, create measurable polarization, even in the absence of 
general relativity. 

5.3.4. Geometrical Optics of Vacuum Waves in a Curved Spacetime 

We first consider the propagation of photons in a vacuum surrounding a black hole. 
These follow orbits called null geodesies, just as material particles follow timelike geodesies. 
The equation of motion can be expressed in a general coordinate system though, in our 
case, Boyer-Lindquist coordinates, by saying that the total derivative of the wave vector 
along the ray vanishes. In index notation, this becomes 

fc"A# = , (5.67) 

cf Eq. ( ^62| ). 



There are essentially three constants of the motion that describe these orbits, an energy, 
— k an angular momentum, k^, and a third quantity known as the Carter constant, Q 
(e.g. Misner, Thorne & Wheeler 1973.) (In fact, we only need the ratios k^/k^.Q/k^ to 
define the rays.) Close to the black hole the rays are strongly curved with the consequence 
that a distant observer, able to resolve a black hole, would be able to see a distorted image 
of the disk behind the hole, apparently hovering above the hole (cf Fig. (2).). The ray 
trajectories are given, in general, by the solution of a set of coupled ordinary differential 
equations that can be partly integrated in terms of elliptic functions (Rauch & Blandford, 
1994). Given a model of the emission, for example of the surface emissivity of a thin 
accretion disk, it is a straightforward, though quite lengthy, exercise to compute the total 
emergent flux and, indeed, the form of the image that would be resolved if the black hole 
could be resolved. 

Now, turn to the propagation of the polarization of vacuum modes in a curved space- 
time, specifically outside the horizon of a Kerr hole (Laor, Netzer & Piran 1990). As 
we have already emphasised, both non-thermal emission (e.g. synchrotron radiation) and 
electron scattering are likely to create polarized sources of radiation. If we just con- 
sider linear polarization for the moment, and the generalization to circular polarization 
is straightforward, then the question that we must answer is "How do we propagate the 
plane of polarization from one point to the next along a curving ray?" . The answer is 
that the electric vector is "parallel-transported" (eg Misner, Thorne & Wheeler 1973). 
What this means is that the unit vector in the direction of the electric field e a changes 
along the ray such that its magnitude and projection onto the direction of the ray remain 
constant. An additional constraint arising from Maxwell's equations in vacuum is that 
the electric vector must be perpendicular to the wave vector. Hence, in index notation, 
we have 

k p e^\ e Q fc Q = 0. (5.68) 

These can be solved consistently to propagate the electric vector along a ray. 
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What is actually done is somewhat different. It turns out that there is another con- 
served quantity, called the Walker-Penrose (1970) tensor, associated with the photon 
spinors (the familiar geometrical object that in this case are associated with light-like 
geodesies, the path taken by photons in vacuum.) This actually involves the electric 
vector and can be used to relate the polarization at the point of emission to that at the 
point of observation directly without having to integrate a differential equation (Connors, 
Stark & Piran 1980). It is then possible to define a transfer function for the polarization 
and to compute the polarized flux given a specific emission model using the propagated 
intensity. (Note that, in propagating the intensity, we must correct for the Doppler and 
gravitational shifts. There is a natural way to do this in general relativity.) 

5.3.5. Geometrical Phase 

The next level of complication is to introduce the plasma into the curved space time. 
Let us do this in two stages. The first stage is to ignore the magnetic field so that the 
local dispersion relation takes the form fl = {u)p{x) + c 2 fc 2 ) 1 / 2 . In this case, the refractive 
index is locally isotropic. This means that the two eigenmodes at a point are degenerate 
and that we have to formulate a rule to connect the polarization from one point along a 
path to the next. 

The wave packets, which travel at the local group velocity, are no longer moving along 

null geodesies, but timelike geodesies instead. It turns out that the Hamiltonian equations 

of motion can be generalized in a covariant manner, provided that one has knowledge of 

the local linearized dispersion relation at every relevant point in space time. Therefore, 

there is a prescription for computing the paths. These can be described by a four velocity 

u a and an acceleration a a = du a /dr with respect to a freely-falling frame, where dr is an 

interval of proper time. The standard relativistic way to handle this is to generalize the 

notion of parallel-transport to Fermi- Walker transport (e.g. Misner, Thorne & Wheeler 

1973), which corrects for the non-null motion of the wave packets. The propagation 

equation becomes 

de a 

— = e p a p u a - u^a a . (5.69) 

or 



Eq. (5.69) reduces to Eq. (5.68) when Eq. (5.67) is satisfied. This provides a natural 
basis in which to discuss polarization propagation and phase changes. 

It is instructive to consider a wave propagating along a twisting optical fiber, with k 
parallel to the local tangent to the fiber. Here again we have gradients in an isotropic 
refractive index. In this case, the unit electric vector, along e, must remain perpendicular 

to the unit wave vector k. When the fiber bends, the change in e must be along k; there 
is no other vector to be involved as the medium is isotropic. Therefore we can write 
down the equation of propagation for the electric vector from first principles. 



This is a limiting case of Eq. (5.6S) 



A good way to visualize what is happening (Berry 1990) is to allow the tangent to 
the fiber to trace out a path on the unit sphere, e is tangent to the sphere and it is 
straightforward to see that rotation angle of e after traversing a complete circuit equals 
the solid angle enclosed by that circuit. If we propagate a linearly polarized wave along 
a twisting fiber, the polarization direction will, in general, be rotated between two points 
where the fiber is parallel. This experiment has been performed successfully (Chiao 
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et al. 1989). (Actually this was under conditions when physical as opposed to geometrical 
optics applies, though the results should be identical.) 

This rotation - essentially a phase change between the two circularly polarized eigen- 
modes - is known as the geometric phase. Geometric phase is a quite general phenomenon 
in physics and analogs are expected to be relevant to wave propagation in a curved space- 
time. The Foucault pendulum provides another example of this general phenomenon. As 
is well known, a Foucault pendulum at latitude £ will only rotate through an angle in 
incrtial space of 2ir(l — sin^), the solid angle traced out by the radial vector on the unit 
sphere, as it is carried around a complete circuit in one day by the spinning Earth. Now 
to see where general relativity may come in, it is helpful to consider a Foucault pendulum 
at the North pole. According to the above discussion, there is no rotation of the plane 
of oscillation according to Newtonian dynamics. However, the tiny dragging of incrtial 
frames effect associated with the spin of the earth leads to an equally tiny rotation of the 
plane of polarization and there have been proposals to measure it. Effects such as these 
would be much larger near a spinning black hole and could also influence the propagation 
of electromagnetic waves. 

5.3.6. General Relativistic Magnetoionic Theory 

The second stage is to reinstate the magnetic field which, on general grounds, is surely 
present. This breaks the degeneracy between the two wave modes. This is akin to 
changing the pivot of a Foucault pendulum from a point attachment to an axle. As far 
as is known, there is no generalization of the Carter constant and the Walker-Penrose 
tensor, to non-null wave modes. The equations for the paths can be integrated and the 
intensity and the polarization can be propagated using the relativistic generalization of 
the approach outlined above (Broderick & Blandford, in preparation.) 

When a linearly-polarized vacuum wave crosses the ergosphere of a rapidly spinning 
hole, there is a contribution to the emergent polarization position angle of order unity 
due to the dragging of inertial frames. We can think of this as a phase difference in 
of order unity in the two circular polarized modes into which the linear mode can be 
decomposed. Now, if we introduce a magnetoactive plasma into the path, then there 
is likely to be a large Faraday rotation per unit length. The total rotation will differ 
by much more than O(l) along different paths that any emitted linear polarization is 
likely to be erased. However, as discussed above, the eigenmodes are not completely 
circular and have an ellipticity 0(XY). What this means is that a systematic phase 
difference O(l) will be introduced between the two modes and that, if the original modes 
arc in phase so that there is no circular polarization, a degree of polarization O(XY) 
will emerge, independent of the reversals of the magnetic field and variations in the total 
Faraday rotation along different lines of sight. 

These matters deserve further attention. 

5.4. Interpretation 

5.4.1. Sgr A* and other Low Luminosity AGN 

X-ray observations of Sgr A* (Baganoff et aL2001) have shown the the luminosity 
is very small (~ 4x 10 33 erg s _1 ) and the spectrum is quite soft. Variability on an 
hour timescalc may also have been seen. (As has been argued elsewhere, Blandford & 
Begelman 1999, this is generally to be expected if most of the mass supplied to the hole is 
blown away in a wind.) This suggests that the density of gas is very low close to the hole 
and opens up the possibility that we may be seeing radio or, more likely, mm emission 
from the ergosphere. Under these conditions, polarization observations could be quite 
diagnostic of the physical conditions. 
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The variable circular polarization discussed above, increasing in degree with frequency, 
has at least three explanations. Firstly, there could be a radius to frequency mapping so 
that the radio photosphere shrinks with frequency and the field gets stronger so that the 
energy of the emitting relativistic electrons also decreases. This leads to an increase in 
the emitted degree of circular polarization (cf §1). Variability studies at high frequency 
should be quite diagnostic. Secondly, the polarization could be due to a flat space 
propagation effect along the lines discussed above. It will be particularly interesting to 
see if the sign of the circular polarization really does not change, as the observations to 
date may suggest. This could, in principle, be related to the angular velocity of the disk 
as outlined above, though quantitatively this seems improbable. 

The third possibility may be the most unlikely, yet it is the most exciting. This is that 
the circular polarization reflect directly the geometry of the ergosphere and be due to a 
general relativistic, propagation effect as outlined above. In the case of Sgr A*, we expect 
the gas at high latitude in the ergosphere to be moving with speed ~ c. We deduce that 
the expected degree of circular polarization is: 

C ~ XY ~ 10 ~ 2 (lo^Fr) (iWc) ( li) 3 ' (5J1) 

It is not impossible that this effect is observable These matter deserve more attention, 
both observational and theoretical, in Sgr A* and other, nearby galaxies. 

5.4.2. Imaging the Ergosphere 

Interesting and timely as these ideas may be, radio observations are probably not likely 
to contribute to a confirmation of the essential features of the Kerr metric until we can 
actually image the ergosphere. At present, as we have remarked, the best resolution has 
been achieved in M87 (~ IOOto). Probably the best prospects lie with sub mm VLBI 
observations of Sgr A*, where interstellar scattering precludes resolving the ergosphere 
at radio wavelengths (Falcke, Melia & Agol 1999). There are also quite futuristic plans 
to develop X-ray interferometry to achieve analogous goals (Cash et <zZ.2000). 

5.5. Summary 

• Black holes are common features of the evolution of massive stars. 

• Massive black holes are commonly found in the nuclei of normal galaxies. Pre- 
sumably they powered active galactic nuclei including quasars and radio sources in the 
past. 

• We have good grounds to be confident in the general theory of relativity and, specif- 
ically, the Kerr metric which describes the curved spacetime around a spinning black 
hole. However, this does not absolve us from the responsibility of testing the theory. 

• We also want to understand how black holes accrete and how they form jets as well 
as their impact on Galactic and extragalactic astronomy. 

• Recent observations suggest that we may be observing radio and mm emission from 
very close to the black hole in Sgr A*. There is consequently interest in developing the 
magnetoionic theory and radiative transfer in a general relativistic environment. 
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